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ABSTRACT 

Wc have obtained Spitzer IRS spectra and MIPS 24, 70, and 160 /im photometry for a volume- 
limited sample of 22 SDSS-selected Low-ionization Broad Absorption Line QSOs (LoBALs) at 0.5 < 
z < 0.6. By comparing their mid-IR spectral properties and far-IR SEDs with those of a control 
sample of 35 non-LoBALs matched in Mj, we investigate the differences between the two populations 
in terms of their infrared emission and star formation activity. Twenty five percent of the LoBALs show 
PAH features and 45% have weak 9.7/im silicate dust emission. Wc model the SEDs and decouple the 
AGN and starburst contributions to the far-infrared luminosity in LoBALs and in non-LoBALs. Their 
median total, starburst, and AGN infrared luminosities arc comparable. Twenty percent (but no more 
than 60%) of the LoBALs and 26% of the non-LoBALs are ultra-luminous infrared galaxies (ULIRGs; 
Lm > 10 12 Z/©). We estimate star formation rates (SFRs) corrected for the AGN contribution to the 
FIR flux and find that LoBALs have comparable levels of star formation activity to non-LoBALs when 
considering the entire samples. However, the SFRs of the IR-luminous LoBALs are 80% higher than 
those of their counterparts in the control sample. The median contribution of star formation to the 
total far-infrared flux in LoBALs and in non-LoBALs is estimated to be 40-50%, in agreement with 
previous results for PG QSOs. Overall, our results show that there is no strong evidence from the mid- 
and far-IR properties that LoBALs are drawn from a different parent population than non-LoBALs. 
Subject headings: galaxies: active - infrared: galaxies — galaxies: evolution — quasars: absorption 
lines — quasars: general 



1. INTRODUCTION 

Supcrmassive black holes are found to be ubiqui- 
tously present at the centers of galaxies with bulges 
and several relationships between the mass of the black 
hole and properties of the spheroid strongly suggest co 



evolution of the two (e.g.. iKormendv fc Richstonel | 1995 
Magorrian et al.lll998t lLaorill998t iGebhardt et al.lboOO ; 



Ferrarese fe Merrittj 12000 1: Kormendv fe Gebhardtl boOl 

H2002T m 



McLure k Dunlod [20021: iTremaine et all 120021) . The 



mechanisms via which the galaxy and black hole reg- 
ulate each other's growth are still unknown. Various 
types of outflows, such as supernova winds and AGN- 
driven winds, have been invoked as the plausible feed- 
back processes responsible for quenching the star forma- 
tion in the host and clearing the gas from the nuclear 
region, and, thus, halting the a ccretion onto the blac k 
hole and limiting its mass (e.g.. lDi Matteo et al.ll2005h . 
However, observational evidence of the e xtent of their in- 
fluence is still s parse and uncertai n (e.g ., lMoe et al.ll2009t 
iBautista etUlkoiOl: iDunn et al.ll2010f) . 

Observations of young, recently fueled QSOs are the 
key to testing this evolutionary model. Particular at- 
tention has been paid to studying the properties of ul- 
traluminous infrared galaxies (ULIRGs; Ljb, > 10 12 £©) 
since they are believed to be powered by both AGN 
and starbursts, although starbursts arc generally be- 
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lieved to be responsi ble for the bulk of the power output 
(I Sanders et al.1 119881 : for review see iSanders fc Mirabel 
I1996D . Theconnection between AGN and ULIRGs is 
suspected by the fact that they are some of the most 
luminous sources in the universe with comparable lumi- 
nosities of Ltoi > 10 12 L Q . In addition, they are both 
associated with s trong infrared emission from dust (e.g., 
IHaas et all 120031 ). The morphologies and dynamics of 
ULIRGs indicat e that these galaxie s are associated with 
galaxy mergers (lArmus et al.lll987t ISanders et al.|[l988t 
Murphy et ail 119961 iVeilleux et al.1 120021: iDasvra et al.l 
20061 ). Similarly, many QSO hosts at z < 0.4 show 



signs of interaction, even some of those that had previ- 
ously been classified as ell i pticals dCanalizo fc Stocktonl 
l200l ICanalizo et al J 120071 iBennert et al.ll2008fl . If star- 
bursts, ULIRGs, and AGN are connected in an evolution- 
ary sequence which was initiated by a galaxy interaction, 
observations of the transition stages of this process are 
necessary to better underst and this connection (see e.g., 
IHonkins et H1l2f)(T7l l2008bl lah. 

BAL QSOs are promising candidates for newly emerg- 
ing optical QSOs. BAL QSOs are a subclass of QSOs 
characterized by broad absorption troughs of UV reso- 
nance lines, blueshifted relative to the QSO's rest frame, 
which ar e indicative of ga s outflows with speeds of up 
to 0.2c tiFoltz et al.1 119831) BA Ls were rigorously de- 
fined by iWevmann et alT(|1991f ) to include only objects 
with broad absorption lines wider than 2000 km s _1 , 
blueshifted past the first 3000 km s" 1 ; however, some 
studies more recently have been more inclusive of the 
wide range of absorption observed and have relaxed that 
criterion to lo wer limits on the ab sorption width of 1000 
km s^ 1 (e.g.. iTrump et all 120061 ). Hydrodynamic mod- 
els show that AGNs are capable of launching such high 
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velocity winds (iMurrav et al.l 119951 iProga et all 120001: 
iGallagher et al.ll2012t ) . Based on the material producing 
BAL troughs, there are at least three subclasses of BAL 
QSOs. The high-ionization BAL QSOs (HiBALs) arc 
identified via the broad absorption from C IV A1549, but 
they might have absorption from other high-ionization 
species such as Lyq, N V A1240, and Si IV A1394 
(jHall et al.H2~002t) . The low-ionization BAL QSOs (LoB- 
ALs), in addition to the lines present in HiBALs, feature 
absorption lines from Mg II AA 2796,2803, Al III, and 
Al II. A very small fraction of LoBALs, called FeLoBALs, 
show absorption in th e rest-frame UV fro m metastable, 
excited states of Fe II ((Hazard et al.lll987t) . 

It is not well understood why only 10%— 30% of the 
optically selected QSOs have BALs dTolea et al.l 120021: 
Hewett fe Folti [20031 : iTrump et all 120061 : iGibson et al.1 
20091). and only about one tenth of these are LoBALs 
(jReichard et al.ll2003af ). Due to the highly obscured na- 
ture and much redder continua of these objects, opti- 
cal identification omits a large fraction of BALs. There- 
fore, although LoBALs are observed in only 1-3% of all 
optically-selected QSOs, they comprise a much higher 
fraction of the QSOs selected at l onger wavelengths 
(jUrrutia et al.ll2?I?)9L IDai et al.l[27)10h . lAllen et al.l pll! ) 
estimate that the intrinsic fraction of BAL QSOs can be 
as high as ~40% when the spectroscopic incompleteness 
and bias against selecting BAL QSOs in the SDSS arc 
taken into account. Hence, BAL QSOs may only be rare 
in optically-selected samples. 

Models attempting to explain their occurrence need 
to account for their obscured nature. Currently there 
are two competing interpretations of the BAL phenom- 
ena: orientation and evolution. On one hand, BAL and 
non-BAL QSOs are thought to derive from the same 
parent populatio n because of the remarkable similar- 
ity in their SEDs (jWevmann et al.lll99 it IGallagher et al.l 
l2007f ). Nonetheless, QSO continua appear to be increas- 
ingly reddened in a sequence going from non-BALs to Hi- 
BAL s to LoBALs (jReichard et alil2003bt iRichards etTaTl 
l2003h . This finding inspired efforts to explain the 
low occurrence of BAL QSOs within the framework 
of the AGN unification model (Antonucci 1993), sug- 
gesting that, due to orientation effects, BALs are seen 
in classical QSOs only when viewed along a narrow 
range of lines of sight passing through the accretion 
disk wind. In this picture, high column density accre- 
tion disk winds of ionized gas are d r iven via resonance 
line a bsorption dMurrav et al.l 119951 : iMurrav fe Chiang! 
[19981: lElvisI [20001 1 This model explains the low oc- 
currence of BALs as a natural consequence of the fact 
that BALs are observed only at a small range of view- 
ing angles. Although BALs are pred ominantly radio 
quiet sources (jStocke et al.l 11984 fl992l ). radio observa- 
tions of the few radio-detected BALs provide a test to 
the orientation of the BAL wind with respect to the ra- 
dio jet. Radio-detecte d BALs are observ e d at a wide 
range of inclinations dBecker et al.l 120001: iGregg et al. 
20001: iBrotherton et all 120061 : iMontenegro-Montes et al. 
20081 : IDiPompeo et alJl2010l ) suggesting that the occur- 
rence of BALs is not a simple orientation effect (e.g., 
IDiPompeo et al.ll20 12f) . Currently it is not clear whether 
or not radio-loud (RL) and radio-quiet (RQ) QSOs arise 
from the same parent population, so it is certainly pos- 
sible that RL and RQ BAL QSOs are different classes. 



An alternative model proposes that BAL QSOs are 
young QSOs caught during a short-lived phase in their 
evolution when powerful QSO-driven wi nds are blow- 
ing away a dusty obsc uring cocoon (e.g ., lHazard et al.1 
H98l iVoit et all Il99l lHall et al.1 [20021 and references 
therein). This model appears to be particularly appli- 
cable to LoBALs since these obj ects are suspected to be 
young or recently r efueled QSOs (jBoroson fe Greenl[l99l : 
iLipari et al.l I1994D and might be exclusi vely associated 
wi th mergers dCanalizo fe Stocktonl [200l. Observations 
by iCanalizo fe Stockton! (|2002h of the only four known 
LoBALs at z < 0.4 at the time showed that: (1) they are 
ULIRGs; (2) they have a small range of far-IR colors, in- 
termediate between those characteristic of ULIRGs and 
QSOs; (3) their host galaxies show signs of strong tidal 
interactions, resulting from major mergers; (4) spectra of 
their hosts show unambiguous interaction-induced star 
formation with post-starburst ages < 100 Myr. Simi- 
larly, studies of FeLoBA Ls, both at low dFarrah et al.l 
[20051 ) and high redshifts (|Farrah et al.ll2007U2010D . sug- 
gest that they are associated with extremely star-forming 
ULIRGs. 

Most rece nt hydrodynam i c simu lations of major galaxy 
mergers by iDeBuhr et all (|2011D show that an AGN- 
driven BAL wind with an initial velocity ~ 10000 km 
s" 1 would lead to a galaxy-scale outflow with velocity 
~ 1000 km s^ 1 , capable of unbinding 10—40% of the 
initial gas of the two merging galaxies. Such AGN feed- 
back could possibly explain the o bserved high-velocity 
outflows in post-st arburst galaxies ([Tremonti et al.ll2007f) 
and ULIRGs (e.g. IChung et alJl2011t IRupke fe VeiUeuxl 
120111: iSturm et al.l 120 111) Further, if the paradigm sug- 
gesting that AGN feedback is responsible for regulating 
the growth of galaxies is correct, LoBALs may be at a 
unique stage where strong outflows are present, yet, star 
formation is still in the process of being quenched. 

Previous studies o f large samples of BA L QSOs ad- 
dressing their SEDs (IGallagher et al.ll2007t) and submil- 
limeter properties (|Willott et al.l 120031) find that BAL 
and non-BAL QSOs are indistinguishable, consistent 
with the model that all QSOs contain BAL winds, and 
their detection depends on viewing angle. However, 
those samples mainly comprise HiBALs and refrain from 
drawing conclusions about LoBALs. Even if the detec- 
tion of BAL troughs in QSO spectra depends on view- 
ings angle, compelling evidence suggests that LoBALs 
are linked to IR-luminous galaxies, with dominant young 
stellar populations and disturbed morphologies. To test 
this possibility, we have undertaken the first multiwave- 
lcngth investigation of a volume-limited sample of LoB- 
ALs. In a series of three papers, we address the na- 
ture of low-redshift LoBALs and their relationship to the 
broader QSO population. In particular, we test the idea 
that LoBALs might be a short, evolutionary stage when 
the AGN has been recently fueled by a merger and the 
ensuing winds are in the process of quenching the star 
formation. 

In this first paper of the series, we present Spitzer IRS 
spectroscopy and MIPS photometry at 24, 70, and 160 
/jm of a volume- limited, statistically-significant sample of 
low-redshift, optically-selected LoBALs. To study their 
star-forming histories, we model the infrared SEDs of 
LoBALs and measure their far-infrared luminosities and 
star formation rates. In the upcoming papers, we will 
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study the detailed morphologies of LoBAL host galaxies 
via HST imaging and the nature of their stellar popula- 
tions via Keck LRIS spectroscopy. 

Our sample of LoBALs and a control sample are de- 
scribed in § 2. Details of the observations are explained 
in § 3. We present the analysis and results in § 4. Dis- 
cussion and summary of results are given in § 5. The 
conclusion is presented in § 6. We assume a flat universe 
cosmology with Ho = 71 km s _1 Mpc -1 , Qm = 0.27, 
and SI a = 0.73. All luminosities in units of the bolomet- 
ric solar luminosity were calculated using Lq = 3.839 x 
10 33 erg s _1 . 

2. SAMPLE SELECTION 
2.1. LoBALs 

We selected a volume-limited sample of 22 LoBALs 
from the Sloan Digital Sky Sur vey fSDSS: lYork et al.1 
[2000]) Third Data Release (DR3; lAbazaiian et al.ll2005D . 
One of the spectral lines which classifies QSOs as LoB- 
ALs is the blue-shifted broad absorption of Mg II A2800. 
Samples of optically-identified LoBALs are only available 
for z > 0.48, when this low-ionization broad absorption 
line is redshifted from rest-frame UV into the spectral 
range of the SDSS. While now there are a vailable SPSS 
catal ogues with LoBALs at z > 0.48 (e.g., iGibson et al.l 
I2009T) . at the time the sample for this project was se- 
lected the most up-to-date ca talog of BAL quasars was 
the one bv I Trump et al.1 ()2006f ). which identified LoBALs 
for z > 0.5. 

In addition to characterizing the FIR emission of LoB- 
ALs, we wanted to study their mid-infrared spectral 
properties, in particular, the PAH emission lines at rest- 
frame 6.2, 7.7, 8.6, 11.3, and 12.8 /j,m and the silicate 
emission at rest-frame 9.7 /xm. Choosing objects at red- 
shifts less than 0.6 ensures that these spectral lines are 
observable within the wavelength range 7—21 /im, allow- 
ing us to use only two of the four channels of the Spitzer 
infrared spectrograph (i.e., SL1 + LL2; see Section [3] for 
details). An upper limit of z < 0.6 also makes it feasible 
to resolve the host galaxy morphologies on HST images 
despite the bright nuclear emission. An attempt to study 
the morphologies of LoBALs at 0.9 < z < 2.0 showed 
that the host galaxies could not be resolved at those red- 
shifts (i.e., Fan; HST PID 10237), which strongly moti- 
vated us to choose a lower redshift upper limit. 

With these considerations in mind, our sample was se- 
lected to include all LoBALs within the r edshift range 0.5 
< z < 0.6 from the ITrump et al.1 (|2006D catalog of BAL 
quasars, drawn from Data Re lease 3 of the Sloan Di gital 
Sky Survey quasa r catal og bv I Schneider et all (|2005l ). In 
the iTrump et all (|2006t ) catalog, QSOs with regions of 
flux at least 10% below the continuum, spanning over a 
velocity range of at least 1000 km s _1 blue- ward of Mg II 
A2800, were identified as LoBALs. They found 457 LoB- 
ALs in the redshift range 0.5 < z < 2.15. Of those, 
only 22 fall within 0.5 < z < 0.6, when we exclude one 
object which is classified as a narrow-line LoBAL and 
one identified as an uncertain FeLoBAL. The volume- 
limited sample of 22 low-redshift LoBALs is listed in Ta- 
ble Q] Note that some of the objects in o ur sample are 
not id entified as LoBALs in the catalog bv IGibson et al.1 
(120091) since the Gibson et al. introduce a new balnic- 
ity index, which is different from the absorption index of 



ITrump et ail (I2006D . The SDSS spectra of the objects 
in our sample, showing the Mg II absorption trough, are 
included in the Appendix. 

2.2. Control sample of non-LoBAL type-1 QSOs 

In order to investigate the relationship between LoB- 
ALs and classical QSOs, we compiled a sample of type-1 
QSOs with sufficient available archival data to be ana- 
lyzed in the same way as the LoBAL sample. We selected 
objects whose spectra did not show Mg II BAL absorp- 
tion. Since we do not have UV spectra to determine 
whether any of these objects are HiBALs, we shall refer 
to them as "non-LoBALs." We drew the control sample 
from t he compilation of quasar SEDs bv lRichards et afl 
(|2006h . who published all available SDSS, Spitzer MIPS, 
and IRAC photometry, as well as x-ray and radio data, 
of the type-1 quasa r s in t he SDSS-DR3 quasar catalog 
by iSchneider et al.l (|2005l ). It is important to match 
the LoBALs to typc-1 QSOs of comparable luminos- 
ity. Limiting the control sample to the same redshift 
range as the LoBALs, i.e., 0.5 < z < 0.6, selected only 
11 non-LoBALs whose absolute Mj magnitudes matched 
only the lower luminosity LoBALs. To avoid the prob- 
lem of having an under-luminous control sample, we 
matched the control sample to the reddening-corrected 
optical luminosity of the LoBALs, not restricting the red- 
shift. Moreover, only fe w LoBALs have been found to 
have radio emissi on (e.g.. iBecker et al.lll997l 120001 120011 : 
iBrotherton et alj|1998t iMenou et al.ll2001l ). which limits 
our selection to the radio-quiet sub-population. 

In summary, we selected all radio-quiet non-LoBALs 
from iRichards et al.1 ()2006l ). which matched the absolute 
Mi magnitude range of our LoBAL sample, —22.4 > 
M t > -25.6 (see Table [TJ. A total of 35 objects within 
0.45 < z < 0.83 satisfied these criteria: 16 from the 
Lockman hole, six from the ELAIS-N1, six from the 
ELAIS-N2 fields in the Spitzer Wide- Area Infrared Ex- 
tragalactic Survey (SWIRE: lLonsdale et aTll2003[ ). and 
seven from the Spitzer observation s of the Extragalac - 
tic First Look Survey area (xFLS: iFraver et al.ll2006l) . 
Tables [TJ [2] and [3] list relevant optical and infrared pho- 
tometry of the comparison sample of non-LoBALs. 

3. OBSERVATIONS AND DATA REDUCTION 

3.1. LoBALs 

3.1.1. Spitzer IRS 

Mid-infrared spectra of 20 of the 22 LoBALs were 
obtained with the Spitzer Infrared Spectrograph (IRS; 
iHouck et al.1 120041) as part of our Cycle 5 GO program 
(Program ID 50792). The IRS observations of the last 
two scheduled objects, SDSS J023102-083141 and SDSS 
J023153— 093333, could not be completed due to the 
depletion of the cryogen and the early commencement 
of the Warm Spitzer mission. All 20 of the objects 
were observed in staring mode with the Short Low first- 
order module (SL1), which covers 7.4—14.5 /im, and with 
the Long Low second-order module (LL2), which covers 
14.0—21.3 /Ltm. Observations in SL1 and LL2, respec- 
tively, consisted of a 6 or 14 s ramp for 2 cycles with two 
nod positions and a 120 or 30 s ramp for one to three 
cycles with two nod positions. The exact ramp dura- 
tions and number of cycles for each object are listed in 
Tablc|4j The slit width was 3'.'7 and 10'.'5 (corresponding 
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TABLE 1 

SDSS DR7 PSF magnitudes corrected for the given here Galactic extinction. 
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in 
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Q 
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01 


16.92 


± 


Q 


00 


16 


92 
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00 


16 


77 


-j- 





1 


16 


89 


-J- 





01 


20 


T1 7001 


8.34-30^4^ 

O- > ,J:.):.}-±.: ) 


g 





577 


-23.40 


0.02 


19.80 


± 





04 


19.45 


± 


() 


01 


19 


26 


-j- 





01 


19 


03 


-j- 





01 


18 


91 


-£ 


() 


05 


21 


T1 70^41 


89-I-3S.3044 


7 





554 


-24.14 


0.04 


19.95 


± 





04 


19.08 


± 





01 


18 


65 


-j- 





01 


18 


23 


-j- 





01 


18 


17 


-£ 


() 


03 


99 


T904^ 33 


9n nni 1 n4 


9 





545 


-24.82 


0.06 


18.43 


± 





02 


18.01 


± 


Q 


1 


27 


79 


-J- 





01 


17 


56 


-j- 





1 


1 7 


64 


_^ 


() 
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Control sample of 


non- 


LoBALs 




































1 


J103236 


.22+580033 


!) 





.687 


-23.20 


0.006 


20.23 


± 


11 


05 


19.85 


± 





.03 


19 


.82 


± 





03 


19 


63 


± 





.03 


19 


.79 


± 





08 


2 


J103333. 


92+582818. 


.8 





.574 


-22.17 


0.007 


20.78 


± 





07 


20.29 


± 





.03 


20 


.40 


1 


0. 


01 


20 


22 


1 





05 


20 


.06 


± 





.14 


3 


J103651. 


94+575950. 


9 





.500 


-23.45 


0.006 


19.48 


+ 





03 


19.05 


± 





02 


18 


79 


1 





02 


18 


60 


1 





02 


18 


.50 


± 





04 


4 


J103721. 


15+590755. 


.7 





.603 


-23.03 


0.008 


20.04 


+ 





.15 


19.43 


± 





.03 


19 


.59 


i 





02 


19 


.49 


-1 





02 


19 


67 


± 





10 


5 


J104210 


25+594253. 


.5 





.675 


-23.17 


0.012 


19.85 


+ 





00 


19.66 


± 





04 


19 


.65 


i 





02 


19 


63 


i 





.04 


19 


.70 


± 





12 


6 


J104526. 


73+595422. 


6 





.646 


-23.85 


0.011 


19.44 


+ 





05 


19.06 


± 





.03 


19 


14 


± 





02 


18 


83 


-1 





.03 


18 


.84 


± 


0. 


05 


7 


J104556. 


84+570747. 








.541 


-23.24 


0.007 


19.21 


+ 





03 


18.95 


± 





02 


19 


.00 


i 





02 


19 


00 


i 





.03 


18 


.99 


± 





04 


8 


J104625. 


02+584839 


1 





.577 


-23.70 


0.010 


19.06 


+ 





02 


18.76 


± 





.03 


18 


.84 


i 





03 


18 


70 


i 





02 


18 


.78 


± 





04 


9 


J104633 


70+571530. 


.4 





.712 


-24.16 


0.008 


19.03 


+ 





02 


18.72 


± 





01 


18 


.77 


i 





02 


18 


75 


i 





01 


18 


.76 


± 





03 


10 


J 104840 


28+563635. 


6 





.714 


-23.86 


0.007 


19.67 


+ 





04 


19.25 


± 





.03 


19 


.20 


i 


0. 


02 


19 


06 


i 





02 


18 


.90 


± 





05 


11 


J104857. 


92+560112. 


3 





.800 


-24.25 


0.009 


19.15 


+ 


0. 


02 


18.81 


± 





04 


18 


82 


i 





01 


18 


95 


■1 





02 


18 


.77 


± 





04 


12 


J105000 


21+581904. 


2 





.833 


-25.53 


0.008 


17.87 


± 





01 


17.67 


± 





.02 


17 


.64 


t 





02 


17 


77 


1. 





.01 


17 


.68 


± 





03 


13 


J105106. 


12+591625 


1 





.768 


-24.78 


0.009 


18.60 


± 





03 


18.22 


± 





.03 


18 


.27 


.1. 





02 


18 


33 


J 





.02 


18 


19 


± 





03 


14 


J105518. 


08+570423. 


.5 





.696 


-24.13 


0.007 


18.87 


± 





02 


18.59 


± 





02 


18 


.71 


i 





01 


18 


73 


d 





02 


18 


68 


± 





03 


15 


J105604 


00+581523. 


.4 





.832 


-24.36 


0.007 


19.02 


± 





03 


18.89 


± 





02 


18 


.81 


i 





03 


18 


93 


i 





02 


18 


.81 


± 





03 


16 


J105959. 


93+574848. 


.1 





.453 


-23.76 


0.006 


18.86 


+ 





03 


18.36 


± 





.03 


18 


25 


i 





02 


18 


05 


i 





02 


17 


.81 


± 


0. 


02 


17 


J160015 


68+552259 


9 





.673 


-24.36 


0.007 


18.90 


± 





03 


18.54 


± 





.02 


18 


.54 


1 





01 


18 


42 


± 





.02 


18 


36 


± 





04 


18 


J160128 


54+544521. 


.3 





.728 


-24.93 


0.010 


18.37 


± 





02 


18.12 


± 





.02 


18 


.06 


± 





02 


18 


05 


± 





.02 


18 


.00 


± 





03 


19 


J160341 


44+541501. 


5 





.580 


-23.23 


0.008 


19.61 


± 





05 


19.26 


± 





.02 


19 


.30 


± 
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19 


19 


± 





.03 


19 


14 


± 





05 


2(1 


J160523 


10+545613. 


3 
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-23.63 
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19.12 
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± 
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.91 
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18 


76 
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.81 
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18 
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.02 


18 
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02 
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18 
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± 





05 


30 


J171334 


02+595028 


3 
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± 





02 


17.88 


± 





.01 


18 


.00 


± 





02 


17 


96 


± 





.02 


18 


.13 


± 





03 


31 


J171736 


90+593011 


4 





.600 


-23.78 


0.021 


19.21 


± 





03 


18.81 


± 





.02 


18 


.90 
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02 


18 


75 
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.02 


18 


.83 
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J171748 


43+594820 


(i 
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The Nature of LoBAL QSOs 



TABLE 2 

2MASS and Spitzer MIPS infrared photometry. 



# 


SDSS Object ID 






J 






H 










hi 


ho 


/l60 






[Vega ma 


g] 


[Vega mag] 


[Vega mag] 


[m,Iy] 


[mjy] 


[mJy] 


LoBAL sample 


1 


J023102. 49-083141. 2 
















15.71 


± 





23 


5.04 ± 0.20 


< 19.63 


< 29.30 


2 


J023153. 63-093333. 5 
















15.92 


± 





25 


1.17 ± 0.17 


< 20.55 


< 22.21 


3 


J025026. 66+000903. 4 


17 


10 


± 





18 






15.42 


± 





15 


10.18 ± 0.20 


86.66 ± 6.68 


86.64 ± 15.02 


4 


J083525. 98+435211. 2 


16 


43 


± 





13 


15.98 ± 


0.20 


15.14 


± 





12 








5 


J085053. 12+445122. 5 


16 


24 


± 





11 


15.98 ± 


0.20 


14.92 


± 





12 


5.68 ± 0.16 


< 17.56 


< 23.17 


6 


J085215. 66+492040. 8 
























1.30 ± 0.16 


< 17.07 


< 21.88 


I 

8 


JOoOOO(.o(-t-400oOU.O 
11(11 Kl O ^/lOO/l O T 

jiui±oi.yo-t-04zy4.£.( 


10 


78 


± 





17 






15.29 


± 





15 


1.90 ± 0.16 
4.52 ± 0.16 


< 16.08 

AO (21 ^ 7 1 a 


< 20.78 

U7 A A -I- 7 OB 
( .44 It ( .yo 


9 


j lUzoUz.oz-foyiyub.D 






























10 


unci no tt i trocn.in o 

J lU51Uz.77+ozoU4y.5 


17 


31 


± 





23 






15.19 


± 





12 


3.28 ± 0.15 


< 18.41 


< 17.80 


1 1 


J 1Uo4U4. 1 o-t-U4zyoy.o 






























12 


T i i o o o o /ii i ,i o o rj n n n 

j i izozz.4i+4ozouy .y 


10 


07 


± 





12 


15.39 ± 


0.14 


14.45 


± 





09 








13 


J il4U4o.Oz-[-ooZ4oy.U 


10 


95 


± 





20 






15.31 


± 





17 








14 


j louyoz.sy-fuiiyou.D 


15 


93 


± 





09 


15.34 ± 


0.11 


14.61 


± 





10 


10.21 ± 0.24 


< 20.93 


< 27.58 


15 


Ti yinnoc; n i on ct n 
J 14UUZO.OO — Ulzyo ( .u 
























1.79 ± 0.20 


< 28.43 


< 30.49 


16 


TI A t n yliC rjC 1 A GQ A 1 A 1 

J 14iy4D.oD + 4Do4z4.o 
























2.02 ± 0.14 


< 17.10 


d/1 7n _(_ "7 "7d 
04. (0 it I . ID 


17 


J 14zo4y .Z4-[-Uozol i . i 


10 


75 


± 





19 


16.03 ± 


0.17 


15.19 


± 





17 


5.85 ± 0.20 


< 20.50 




18 


J 14zyz ; .zo + ozoo4y.o 


10 


27 


± 





10 


15.59 ± 


0.12 


14.85 


± 





09 


8.84 ± 0.17 


<C 14.17 


< 25.28 


19 


T1C1 /IOK 1 "7 1 Q7C11 fl T 


10 


10 


± 





08 


15.44 ± 


0.12 


14.38 


± 





08 


20.06 ± 0.21 


1 1(1 _(_ CCd 

110. 2D ± D.DO 




20 


TiTnnin OO i gncc/ic o 
J 1 i / UUlU.oo + oyoD4o.o 
















15.65 


± 





21 


4.91 ± 0.16 


40.32 ± D.4o 




21 


J 17Uo41.ts2+ oooy44.7 


10 


82 


± 





17 


16.02 ± 


0.20 


15.44 


± 





17 


8.83 ± 0.17 


< 17.54 




22 


J 204ooo.20 — 001104.2 


16 


84 


± 





15 


15.80 ± 


0.15 


15.24 


± 





15 


5.31 ± 0.18 


< 15.23 


















Control sample of 


non 


-LoBALs 






1 


103236.22+580033.9 


18.46 0.26 


17.59 0.28 




0.00 




1.16 ± 0.02 


< 4.41 ± 1.47 


< 12.46 ± 4.15 


2 


103333.92+582818.8 
















16.57 


± 





15 


0.74 ± 0.02 


< 4.52 ± 1.51 


< 13.12 ± 4.37 


3 


103651.94+575950.9 


16 


96 


± 





19 


16.05 ± 


0.20 


15.40 


± 





18 


4.59 ± 0.02 


27.97 ± 2.14 


< 12.19 ± 4.06 


4 


103721.15+590755.7 


18 


70 


± 





2:5 






16.90 


± 





20 


2.13 ± 0.02 


< 4.36 ± 1.45 


< 14.41 ± 4.80 


5 


104210.25+594253.5 


18 


45 


± 





24 














1.08 ± 0.02 


< 3.92 ± 1.31 


< 17.06 ± 5.69 


6 


104526.73+595422.6 


17 


78 


± 





17 


16.64 ± 


0.14 


16.11 


± 





15 


4.38 ± 0.02 


8.20 ± 2.53 


< 12.48 ± 4.16 


7 


104556.84+570747.0 


17 


80 


± 





14 


16.94 ± 


0.18 


16.28 


± 





15 


2.72 ± 0.02 


< 4.46 ± 1.49 


< 11.76 ± 3.92 


8 


104625.02+584839.1 


17 


97 


± 





17 


16.87 ± 


0.19 


16.05 







13 


3.88 ± 0.02 


< 3.35 ± 1.12 


< 7.79 ± 2.60 


9 


104633.70+571530.4 


17 


57 


± 





12 






16.10 


± 





13 


1.84 ± 0.02 


< 4.75 ± 1.58 


< 9.52 ± 3.17 


10 


104840.28+563635.6 


17 


80 


± 





14 


16.85 ± 


0.15 


16.42 


± 





19 


2.66 ± 0.02 


< 5.78 ± 1.93 


< 11.93 ± 3.98 


11 


104857.92+560112.3 


17 


55 


± 





13 






16.54 


± 





23 


1.37 ± 0.02 


27.09 ± 2.05 


< 12.11 ± 4.04 


12 


105000.21+581904.2 


16 


80 


± 





16 


16.20 ± 


0.24 


15.46 


± 





18 


4.62 ± 0.02 


< 4.62 ± 1.54 


< 12.93 ± 4.31 


1.'! 


105106.12+591625.1 


17 


17 


± 





11 


16.62 ± 


0.13 


15.34 


± 





15 


5.39 ± 0.02 


23.33 ± 2.15 


< 10.41 ± 3.47 


14 


105518.08+570423.5 


17 


01 


± 





10 






16.36 


± 





17 


3.55 ± 0.02 


< 5.49 ± 1.83 


< 9.51 ± 3.17 


15 


105604.00+581523.4 


17 


28 


± 





13 


16.87 ± 


0.21 


16.23 


± 





22 


4.09 ± 0.02 


< 3.70 ± 1.23 


< 11.50 ± 3.83 


10 


105959.93+574848.1 


10 


05 


± 





10 


16.08 ± 


0.22 


15.48 


± 





10 


9.04 ± 0.03 


17.31 ± 2.45 


< 11.23 ± 3.74 


17 


160015.68+552259.9 
















15.70 


± 





22 


4.48 ± 0.02 


17.3733 ± 2.54 


< 11.66 ± 3.89 


IS 


160128.54+544521.3 


10 


78 


± 





24 






16.13 


± 





.30 


12.93 ± 0.02 


38.14 ± 2.57 


< 11.68 ± 3.89 


19 


160341.44+541501.5 


17 


22 


± 





83 














1.18 ± 0.02 


< 6.78 ± 2.26 


< 11.01 ± 3.67 


20 


160523.10+545613.3 
















16.07 


± 





27 


4.74 ± 0.02 


7.88 ± 1.81 


< 10.78 ± 3.59 


21 


160630.60+542007.5 
























5.56 ± 0.02 


< 4.77 ± 1.59 


< 10.17 ± 3.39 


22 


160908.95+533153.2 












16.22 ± 


0.20 










3.03 ± 0.02 


< 4.92 ± 1.64 


< 11.28 ± 3.76 


23 


163031.46+410145.6 
















15.69 


± 





24 


2.17 ± 0.02 


12.03 ± 1.69 


< 4.47 ± 1.49 


24 


163135.46+405756.4 


17 


33 


± 





82 


16.21 ± 


0.24 










4.3 ± 0.02 


10.98 ± 1.99 


< 4.57 ± 1.52 


25 


163143.76+404735.6 


17 


13 


± 





28 






15.96 


± 





27 


3.9 ± 0.02 


13.91 ± 1.88 


< 4.71 ± 1.57 


20 


163352.34+402115.5 


17 


10 


± 





27 














2.86 ± 0.02 


13.88 ± 2.23 


< 4.56 ± 1.52 


27 


163502.80+412952.9 


17 


08 


± 





22 


16.30 ± 


0.21 


15.22 


± 





10 


3.92 ± 0.02 


< 5.13 ± 1.71 


< 5.30 ± 1.77 


28 


163854.62+415419.5 


17 


41 


± 





33 














2.91 ± 0.02 


< 4.21 ± 1.40 


< 4.54 ± 1.51 


29 


171126.94+585544.2 
























3.45 ± 0.07 


29.10 ± 2.24 


98.32 ± 11.27 


30 


171334.02+595028.3 


17 


12 


± 





27 


16.50 ± 


0.24 


15.47 


± 





20 


5.38 ± 0.07 


< 1.99 ± 0.66 


< 20.21 ± 6.74 


31 


171736.90+593011.4 


17 


63 


± 





35 














6.38 ± 0.05 


19.81 ± 1.46 


< 12.33 ± 4.11 


32 


171748.43+594820.6 


17 


20 


± 





33 














3.04 ± 0.04 


10.34 ± 1.30 


< 7.72 ± 2.57 


33 


171818.14+584905.2 
























4.06 ± 0.06 


30.43 ± 2.26 


< 7.05 ± 2.35 


84 


172104.75+592451.4 


17 


39 


± 





8,2 






16.19 


± 





29 


4.44 ± 0.06 


12.56 ± 2.10 


< 18.29 ± 6.10 


35 


172414.05+593644.0 
















15.75 


± 





24 


1.81 ± 0.06 


< 3.58 ± 1.19 


< 16.80 ± 5.60 
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TABLE 3 

Control sample of non-LoBALs: MIPS data fields and archival IRAC photometry. 



# 


SDSS Object ID 




MIPS field 




G"Jy) 


(jjjy) 


/S.8 
(/iJy) 


J8.0 
(fj,3y) 


4 
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Q 
■ 1 
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2 


103333 


,92+582818. 
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Q 

o 
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1 07Q 7 4-110 
j-U 1 y . 1 zn 1 J- .u 
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4 
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7 
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5 
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5 
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u 
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u 
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1 1 Q7 9 4- c: q 
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7 
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n 

u 
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± 2.5 


7 4- O Q 
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009 7 -1- 7 


o 
o 
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I 
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■ > 
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4 


SWIRE-Lockman 
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U 
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14 
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D 
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15 
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4 
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18 
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3 
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19 
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5 
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20 
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3 
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±101 .1 x y . u 


21 
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D 
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29 
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2 
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23 
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Q 
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24 


163135 
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4 
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582.9 
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ft 
u 


SWIRE-ELAIS N2 
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1 001 4- Q 4 
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26 


163352 
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5 
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460.8 
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578.2 ± 3.4 


797.3 ± 9.7 


1083.8 ± 7.3 


27 


163502 


,80+412952 


9 


SWIRE-ELAIS N2 


901.4 


± 3.2 


1186.5 ± 4.9 


1521.2 ± 10.6 


1961.0 ± 9.2 


28 


163854 


,62+415419, 


5 


SWIRE-ELAIS N2 


532.0 


± 2.3 


673.3 ± 3.2 


850.5 ± 7.9 


1131.0 ± 6.7 


20 


171126 


,94+585544, 


2 


xFLS 




641.4 ± 65.1 


813.2 ± 82.4 


951.0 ± 101.2 


1216.0 ± 124.6 


30 


171334 


,02+595028, 


3 


xFLS 




763.2 ± 76.8 


1015.4 ± 102.3 


1456.1 ± 148.7 


1732.3 ± 175.2 


31 


171736 


,90+593011, 


4 


xFLS 




457.5 ± 46.7 


613.3 ± 62.0 


771.2 ± 84.2 


1366.4 ± 138.8 


32 


171748 


,43+594820, 


6 


xFLS 




647.8 ± 65.3 


898.1 ± 90.9 


1198.5 ± 124.3 


1570.8 ± 160.1 


33 


171818 


14+584905 


2 


xFLS 




675.0 ± 68.3 


869.1 ± 88.1 


1152.5 ± 122.7 


1493.2 ± 152.2 


34 
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,75+592451 


1 


xFLS 




362.1 ± 36.9 


453.4 ± 46.6 


642.1 ± 70.0 


853.8 ± 89.2 


35 


172414 


,05+593644 





xFLS 




591.2 ± 59.5 


751.8 ± 76.0 


853.0 ± 89.0 


1124.8 ± 114.8 



to 24 kpc and 67 kpc at 2=0.55), for the SL1 and LL2 
orders, respectively. 

We used the pipeline coadded, non-subtracted post- 
BCD frames. The bad pixels in the images were removed 
with the interactive IDL procedure IRSCLEAN1.9 (ver- 
sion 1.7). Sky subtraction was achieved by subtracting 
the two-dimensional image at one nod position from the 
other nod position of the corresponding order. One- 
dimensional spectra were extracted with the Java-based 
Spitzer IRS Custom Extraction software (SPICE; ver- 
sion 2.3 Final) using the default parameters of the opti- 
mal extraction method for point sources. We used IDL to 
combine the two nod position in each order and then to 
combine the spectra from the two orders into a continu- 
ous spectrum spanning the wavelength range of 7.4—21.3 
[j,m: first, the spectrum at each nod position was median- 
smoothed and interpolated on a uniform wavelength grid; 
then, the two nod positions for each order were averaged; 
and finally, the two orders (e.g., SL1 and LL2) were con- 
catenated by averaging the spectra in the overlapping 
region. The combined mid-IR spectra of the 20 LoBALs 
observed with the IRS are shown in Fig. [TJ 

3.1.2. Spitzer MIPS 

Far infrared photometry at 24, 70, and 160 (im was 
obtained with the Multiband Imaging Photometer for 



Spitzer (MIPS; iRieke et al.l 12001 as part of our GO 
program (ID 50792). Although the entire sample of 22 
LoBALs was scheduled for MIPS observations, only 17 
targets could be observed prior to the cryogen depiction 
and the commencement of the Warm Spitzer mission. 
The 24 and 70 /im observations were obtained using the 
small-field default resolution photometry mode, while for 
the 160 /Win observations we used the small-field enhanced 
resolution mode. Typical observing modes were one 3 s 
cycle at 24 /zm (48 s integration time), one 10 s cycle 
at 70 /im (126 s integration time), and four 10 s cy- 
cles at 160 /mi (179 s integration time). The standard 
tasks of cosmic-ray removal, imagc-coaddition, and dark 
subtraction were carried out by the automated MIPS 
data reduction pipeline at the Spitzer Science Center. 
Our data reduction started from the pipeline BCD files, 
which were assessed to be of sufficient quality. For the 
24 /im observations, final mosaic images were constructed 
with the M Osaicker and Point source E xtractor software 
(MOPEX: iMakovoz fc Marleaul [2005h after flat-fielding 
and background correction. At 70 and 160 /im, we used 
the pipeline-filtered BCDs to construct mosaics for most 
of the sources. Aperture photometry at 24, 70, and 160 
/im was performed with IDL, using an aperture radius 
of 13", 35", and 48", respectively. The MIPS fluxes are 
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TABLE 4 

LoBALs: Sptizer MIPS and IRS observing log. 



4t 
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3 
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1 


29 


30 


3 
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1 
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6 


1 


13 


14 


1 


29 


5 
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14 


1 


29 


30 


3 
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6 
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14 


1 


29 
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1 


244 


7 


J085357.87+463350.6 
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14 


1 


29 


30 


3 
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8 
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6 


2 
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30 


2 
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10 
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listed in Table [5J All of the 17 sources observed with 
MIPS were detected at 24 microns, but only four were 
detected at 70 /im and only three at 160 /im. Small 
field photometry with MIPS 160 ^m, in particular, has 
the problem that the filtering step of the data reduc- 
tion leads to a loss of flux (Sajina 2012, in preparation) . 
The m anual states this loss is 10% (see also lSaiina et al.l 
21 )()>;. but it can be as high as 30—50%. Three sigma up- 
per limits at 70 fxxa and 160 /im were estimated from the 
standard deviation images associated with the mosaics. 
All fluxes have been corrected for the finite aperture size 
by multiplying by a correction factor of 1.16 at 24 /im, 
1.22 at 70 fim, and 1.601 at 160 /im (aperture corrections 
from the MIPS Data handbook). The systematic uncer- 
tainties are 4% at 24 /mi, 5% at 70 /zm, and 12% at 160 
fim. The uncertainties listed in Table [2] for the LoBAL 
sample are from the aperture photometry calculation. 

3.1.3. Near- Infrared and Optical photometry: SDSS and 

2MASS 

All of the 22 LoBALs in the sample have available 
SDSS Data Release Se ven (DR7) and Two Micron All 
Sky Survey (2MASS; ISkrutskie et "all l2006h photome- 
try pu blished in the quasar catalog of iSchneider et all 
([20101 . Table [1] lists the best SPSS PSF ugriz AB mag- 
nitudes (ugriz: lFukugita et al.lll996lh corrected fo r Galac - 
tic extinction using the map of ISchleeel eTaTl (fl998lh 
We impos e more conservat i ve lim its on the near-infrared 
data than Schneid er et al.l (|2010f) by considered reliable 
only 2MASS magnitudes with Photometry Quality Flag 
(ph_qual) A, B, or C and Read Flag (rd_flg) 1, 2 or 3, 
which ensure measurements with signal-to-noise greater 
than five and measurement uncertainty le ss than 0.2. 
The m eaning of the flags can be found in iCutri et al.l 
(|200l . Table [2] lists the 2MASS J (1.25 /xm), H (1.65 
/urn), and K s (2.16 /im) Vega magnitudes of the LoBALs. 



3.2. Control sample 

3.2.1. SDSS, 2MASS, Spitzer IRAC and MIPS 
photometry 

We use the SDSS DR7 u, g, r, i, z and 2MASS J , 
H, K s photometry published by ISchneider et al.l (|2010f ) 
together with the IRAC 3.6, 4.5 , 5.8, 8.0 ixm, and the 
MIPS 24 pirn data published bv iRichards et al.l (|2006f) . 
All the Spitzer data are taken from the xFLS and 
SWIRE ELAIS-N1, ELAIS-N2, and Lockman Hole 
areas. To obtain 70 and 160 /jm fluxes, we performed 
aperture photometr y with IDL on the pr ocessed mosaics 
provided online by iFraver et al.l (120061 ) for the xFLS 
field 4 and bv lLonsdale eFall ((20031 ) for the three SWIRE 
fields 5 . The 3a upper limits at 70 and 160 were 
estimated from the uncertainty images. Table [3] lists the 
fields from which the MIPS photometry for individual 
objects was extracted. 



4. ANALYSIS AND RESULTS 

4.1. MIR spectral features of LoBALs 

The combined SL1 and LL2 orders covers an observed 
wavelength range X b s = 7.4 — 21.3 fim, which at rcdshift 
0.5—0.6 translates to rest-frame A res t « 5 — 13 /im, after 
removing the noisy region on the red side. LoBALs ex- 
hibit a wide range of mid-IR spectral properties (Fig.[T]). 
About one third (40%) of the spectra are featureless low- 
signal-to-noise continua, nearly half (45%) of the objects 
show silicate emission near 10 /zm and a quarter (25%) 
have at least one of the PAH emission lines, with con- 

4 xFLS data from 



<http: / /data. spitzer. caltcch.edu/popular/fls/cxtragalactic_FLS/enhanccd_MIP! 

5 SWIRE data from 
<http: / /swirc. ipac.caltcch.edu/swirc/astronomers/data_access.html> 



8 



Lazarova et al. 



(f) 

CD 
Q 

X 
7J 



15 r(3) 025026 
1 r 
5 r 



(2') -70341 



(17) 142649 



(14) 130952 



(18) 142927 



[-(19) 161425 



(5) 085053 



(12) 1 12822 




(6) 085215 



S 97 =0.81 



S g7 =0.75 - 
PAH - 



S g7 =0.64 - 

PAH : 



S 9 7 =0.56 j 



S 97 =0.43 J 



PAH J 



S g7 =0.34 



S o7 = 0.34 - 



PAH 



Si 9.7 




iu r 1 

Rest wavelength [/zmj 




io r 1 
Rest wavelength [/xm J 



Fig. 1. — Low-resolution IRS spectra of the LoBALs, plotted as flux density in units of mjy vs. rest-frame wavelength in /im. The 
spectra are in solid black, with la errors in gray. Plotted with a long-dash green line is the spline-interpolated continuum used to estimate 
the silicate strength. Vertical dotted lines at 6.2, 7.7, 8.6, 11.3, and 12.8 (im indicate the fiducial peaks of prominent PAH features. The 
objects in the left column show silicate emission at 9.7 /an, listed in order of decreasing strength from top to bottom. Note that 114043 
and 085053 show apparent redshifted silicate emission, the strength of which could not be measure because the line is truncated and MIPS 
photometry was not available to constrain the SED in the FIR. Also note that 085215 shows possible silicate absorption. The top three 
objects in the right columns show the strongest PAH features. All objects labeled 'PAH' in the top right corner show at least one of the 
PAH lines. The displayed spectra have been median smoothed with a boxcar of five. 
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current PAH and silicate emission present in one tenth 
(10%) of the sample. 

4.1.1. PAH emission features 

The emission from polycyclic aromatic hydrocarbon 
molecules (PAHs), producing prominent lines in the 
mid- IR peaking at 3.3 , 6.2, 7.7, 8.6 11.3 and 12.8 
Aim (jGillett et all 119731: iLeger fc Pugetl 119841) . is pow- 
ered by moderate UV radiation. PAHs are observed in 
photo-dissociation regions (PDRs) where young bright 
stars are contiguous with dense molecular clouds, condi- 
tions found in star- forming re g ions and reflection neb - 
ulae (e.g., iDulev et al.l 119911: iVerstraete et al.l I1996D . 
The correlation between PAH emission an d star forma- 
tion rates in star-fo rming galaxies (e.g., IGenzel et al.l 
[19981 iRoussel et all [20001 iDale fe Heloul |2002j) estab- 
lished PAHs as tracers of star formation (but see 
I Haas et aLl (|2002D for a counter argument) . On one hand, 
there is evidence that the PAH molecule s are destroyed 
by the extreme UV and X-rays radia t ion dPuget fe Legerl 
19891 IVoitl IT9921 lAllain etldl [19961 IGenzel et al.l fl998F 



Smith et al.H2007t 10 'Dowd et alJI2009l:IHunt et aLll2010F) . 

so PAH-derived star formation may underestimate the 
actual activity. In AGN, PAH features are absent (e.g., 
Roche et al.|[l99lt IE. Le Floc'h et al.l[200l . weak (e.g., 



Laurent et al.ll20 00). or have low equivalent widths (e.g., 
Clavel et al.l |2000T T implying destruction or inability of 
the nuclear radiation to excite the aromatic feature. On 
the other hand, the AGN may enhance the PAH emission 
since the nuclear continuum contributes ample flux in the 
UV, which may modify the grain d istribution and d i- 
rectly excite the PAH emission (e.g., ISmith et al.l[2007h . 
Thus, adopting PAH emission as a star formation tracer 
in AGN hosts should be done with caution. 

Studies of the spatial distribution of the PAH emission 
in nearby AGN find that the aromatic em ission arises 
in an extended circum-nuclcar region fe. g.. ICutri et aT 



1984 . < 1 kpc) or in the galactic disk (lLaurent et al 



2000f ). which toget her with the low equivalent widths o f 



PAH features (e.g.. lRoche et al.lll991llciavel et al.ll2000t l 
implies that PAHs in AGN are predom inantly excited 
by star-formation (e.g., IShi et all I2007L and references 
therei n) . This is also supported by ISchweitzer et al.l 
(|2006| ) who find the same ratio between the 7.7 /im PAH 
and the 60 /zm luminosity in a sample of 27 PG QSOs 
and in starburst-dominated ULIRGs, suggesting that the 
starburst is producing all of the QSO FIR emission. 

In our sample of 20 LoBALs with available IRS spec- 
tra, we detect weak PAH features only in five of the 
objects. The PAH emission is strongest for the three 
objects shown at the top of the right column in Fig. [TJ 
The entire complex of PAH lines is seen only in SDSS 
J170010+395545, while PAHs arc concurrent with sili- 
cate emission only in two sources. We confirm the previ- 
ously observed low incidence of PAH emission in AGN, 
in general, to be true also for LoBALs, in particular. 

The average MIR spectrum of all LoBALs (Fig. EJ}, ob- 
tained by averaging the signal-to-noise- weighted individ- 
ual spectra after normalization to the 6 /im continuum 
flux, shows very weak PAH emission at 6.2, 11.3 and 12.8 
/xm. Plotted in the same figure are the average spectra of 
LoBALs grouped according to shared MIR spectral char- 
acteristics, i.e., objects with PAH emission, those with 
silicate emission, and those showing neither silicate nor 



PAHs. PAH emission at 6.2 /xm becomes more promi- 
nent in the average spectrum of the LoBALs with PAHs. 
In the average spectrum of LoBALs which did not oth- 
erwise show individual PAH features, we detect the 12.8 
/im line. We note that the feature at 12.8 /im might be 
a blend of the PAH at 12.7 /im and the low-excitation 
fine-structure emission line of [Ne II] 12.8 /im, which is 
dominant in H II regions and used as SF tracer. In QSOs, 
however, [Ne II] could also aris e in th e narrow line region 
of the AGN and lVeilleux et al.l ((2009ft show that the star- 
burst contributes at most 50% to its flux. In addition, 
we subdivided the 15 LoBALs that had both IRS spectra 
and FIR constraints into those with AGN infrared lumi- 
nosity greater than and less than their starburst lumi- 
nosity (see section l4~3l for details on the measurement of 
these quantities) . The two averages are shown in Fig. [2j 
denoted Lff < Lfg N and iff > Lfg N , respectively. 
We note that in the cases of FIR upper limits for objects 
with Ljp > Lf^ N , it is not certain whether Lj^ is less 
than or greater than Lf^ N . 

The PAH lines are stronger in the composite of ob- 
jects with iff > Lfg N than in the "PAH average," 
suggesting that the FIR emission, indeed, arises from 
star formation rather than from a very extended torus, 
for example. The average of all LoBALs is remarkably 
similar to the subgroups of those without PAHs and 
the one with silicate emission. For comparison, we also 
plot an average sp ectrum of the five type-1 QSOs in the 
iHiner et alj (|2009T ) sample which fall within the redshift 
range 0.5 < z < 0.63. The composite of all LoBALs is 
similar to that of their type-1 QSOs. 




6 a 10 

Rest wavelength [/im] 

Fig. 2. — Signal- to- noise- weighted average spectra of the LoB- 
ALs in our sample, grouped according to shared characteristics. 
The average spectrum of all LoBALs (solid black line) was ob- 
tained by averaging the individual spectra after normalizing their 
flux density by the average flux density from 5—6 fim. With a 
solid orange line is the average spectrum of the objects which show 
no or very weak PAH features. In solid red is the average of the 
LoBALs with PAH features. In solid green we plot the average of 
the LoBALs with silicate emission at 10 /jm. In solid purple is the 
average of the LoBALs with starburst luminosity higher than the 
AGN luminosity from 8—1000 (im, while in solid blue with star- 
burst luminosity less than the AGN luminosity. For comparison, 
in solid gray at the bottom we also plot the a verage spectrum of 
the five type-1 QSOs in the Hincr et al. (2009) sample which fall 
within the redshift range 0.5 < z < 0.63. Vertical dotted lines at 
6.2, 7.7, 8.6, 11.3, and 12.8 (im indicate the position of prominent 
PAH features. The PAH lines are stronger in the composite of ob- 
jects with Lj^ > Lf^ N than in the "PAH average," suggesting 
that the FIR emission, indeed, arises from star formation rather 
than from a very extended torus, for example. 



10 



Lazarova ct al. 



4.1.2. Silicate emission feature at 10 [im 

Silicate compounds comprise the majority of the inter- 
stellar dust and manifest themselves in the MIR via two 
main features centered around 10 /jm and 18 /xm. The 
wavelength coverage of our IRS spectra allows us to in- 
vestigate the silicates peaking at 9.7 iim which arise from 
the s tretching of the Si— O bond (|Knacke fc Thomson! 
I1973D - Due to low dust temperatures of order T ~ 100 
K, this feature is also seen t o peak at longer wavel engths, 
up to Xsi ~ 11 /um (e.g. , IZakamska et all 120081 ) . The 
exact peak of the feature is speculated to be sensitive 
to grain size and composition, with larger dust grains 
and crystalline dust (as opposed to amorphous dust) 
causing the feature to peak at lon ger wavelength s (e.g. , 
IBouwman etai] 120011 ). However, INikutta et ail (|2009f ) 
show that the flat-topped silicate emission peaks ob- 
served in several type-1 sources can be explained by sim- 
ple radiative transfer effects in standard Galactic dust in 
a clumpy torus model. 

Most studies of silicate detections in AGN show that 
there is a notable dependence of the silicate strength on 
the optical classification of the AGN. In the framework of 
the unification model, silicate emission arising from the 
dust torus surrou nding the accr e tion d isk was predicted 
by the models of iPier fc Krolikl (Il993f ). If we attribute 
the appearance of the silicate features in AGN to view- 
ing angle, then type-1 AGN are expected to show silicate 
in emission and type-2 objects in absorption. In fact, 
most observations of AGN in which the feature is de- 
tected support this interpretation. Type 2 quasars (e.g. , 
IZakamska et alj2008tlLacv et al. 2007; H iner et al.l2009t) 
and Scyfert 2 galaxies (|Hao et al.l I2007T ) are almost ex- 
clusively characterized by silicate absorption, while both 
weak silicat e emission and absorption are found in type-1 
QSOs (e.g.jSiebenmorgen et all 120051: lHao et al.ll2005bfc 
Schweitzer et al.l 120081: iHiner et al.l 120091: lLandt et al.l 



2010D and Seyfert 1 galaxies (jHao et al.ll2007l ) 



However, silicate emission or absorption is not ubiqui 
tously present in AGN. The occasio nal detection of sil 

e^JT 

ooa 



icate emission in type-2 AGN 
Hao et all l2005b|: iTeplitz et all l200e 



ISturm et al. 
lLacv et al.l 



2005 



2007 



.Hiner et al.TT2009f) and silicate absorption in ty pe-1 ob 
iects fe.g.. lWeedman et al.ll2005l : lHao et aT1l2007l) , as well 
as the absence of the silicate feature in many AGN, chal- 
lenges the universality of the orientation model and sug- 
gests a more complicated scenario. For instance, with the 
exception of two objects, one mini-LoBAL a nd one red 
QSO, all of the type-1 and type-2 QSOs in the lHao et al.l 
(|2007f ) sample exclusively show weak silicatein emission 
which argues against orientation-dependent silicate emis- 
sion i n QSOs, at least. Models of clumpy torus geome- 
tries (jNenkova et all 120021 l2008a|TB) an d /or larger dust 
grain sizes (e.g., lLaor fc Drainel 119931 : iMaiolino et al.l 
l2001t ) offer a solution to these discr epancies, but can onl y 
reproduce weak silicate features. INikutta et al.l (|2009f ). 
for instance, show that clumpy dust geometry of the ob- 
scuring region can explain both why the 10 /im feature 
is not seen in deep absorption in any AGN and why it 
has been detected in emission in type-2 sources. 

In order to compare our LoBALs to other studies, we 
calculate the silicate strength at 9.7 //ni, $9.7, as: 



'9.7 



In 



where / pea fc(9.7//TO) is the observed flux density at the 
peak of the silicate feature and / ctm t(9.7^m) is the con- 
tinuum flux density interpolated below the peak of the 
emission line. The common challenge of determining the 
underlying continuum is further exacerbated by the lim- 
ited wavelength range of our IRS spectra (i.e., rest- frame 
5.0—14.0 fim). As is apparent in Fig. [TJ some of the sil- 
icate features have a truncated red wing and are often 
rcdshifted. Fortunately, the MIPS photometry allows us 
to extend the range of the IRS spectra by modeling the 
SEDs, which includes the silicate emission amplitude as a 
free parameter in the fit (see Section l4T2"j) . Using the over- 
all SED model, we d etermine the continu um by following 
the fitting recipes of iSpoon et al.l (|2007| ) for continuum- 
dominated sources and interpolate the local mid-IR con- 
tinuum over the range 5.0—31.5 /xm by fitting a cubic 
spline to the 5.0-7.5, 14.0-14.2, and 26.1-31.5 yum con- 
tinuum regions. 

In our sample of LoBALs, silicate is present exclusively 
in emission and is detected in nine of the 20 objects, rang- 
ing from S9.7 = 0.34—0.81. Hence, LoBALs exhibit the 
weak silicate emission typ i cal of other type-1 QSOs (e.g., 
Siebenmorgen et al.l 120051 lHao et al.l | 2005bt lHaas et ahl 
20051 : iClearv et al.l 120071 : Miner et all 120091 ). However, 
we note that the detection/non-detection of silicate dust 
emission might be a function of the signal-to noise ra- 
tio of the data (see Section I4.5p . The strength of the 
feature is individually noted in the bottom right corner 
of each spectrum in Fig. [T] The interpolated contin- 
uum is over-plotted with a dashed green line. Although 
apparently present, the strength of the feature could 
not be estimated in two objects, SDSS J112822+482308 
and SDSS J114043+532439, due to lack of MIPS pho- 
tometry at long wavelengths and, consequently, an SED 
model which allows better determination of the contin- 
uum. We note that SDSS J085215+492040 shows pos- 
sible broad, weak absorption, but the IRS spectrum for 
that object has very low signal-to-noisc ratio and our 
inability to measure the feature precludes us from fur- 
ther speculation. Silicate absorption features is never 
observed to peak at wavelengths longer than 9.8 /im (e.g., 
INikutta et al.ll200"9l ). while the peak of this apparent dip 
is located at ~11 /im. 

4.2. SED modeling 

In general, the infrared emission in galaxies hosting 
an AGN is a combination of (1) starburst emission 
from dust heated by the U V flux of O and B stars 
in active star-forming regio ns ([Devereux fc YoundlT990l: 
IDevereux fc Hameedlll997T) . (2) the AGN emission from 
the dusty torus reprocessing the accretio n disk contin- 
uum and re-radiating it in the infrared (jPier fc Krolikl 
Il993h . (3) diffuse, ambient cold dust (cirrus) emis- 
sion illuminated by the interstellar radiation field (e.g., 
iRowan-Robinson fc Crawfordl I1989D . and (4) infrared 
emission from evolved stellar populations (jKnapp et al.l 
[19921: iMazIei fc de Zottilfl99l . 

We fit the optical-through-FIR SEDs of the LoBALs 
and the control sample of non-LoBALs with a multi- 
component empirical model, which allows us to disentan- 
gle the two major sources that power the FIR emission, 
i.e., starburst and AGN. In our phcnomcnological SED 
modeling approach, we assume that cirrus and evolved 
stellar populations have negligible contribution to the 
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FIR power budget, and account only for the starburst 
and AGN contributions. 

The available optical photometry of our LoBAL sample 
is listed in Table [TJ the MIPS photometry in Table El and 
the IRS spectroscopy is plotted in Fig. [TJ The photom- 
etry of the control sample is given in Tables [TJ , El and EH 
The f itting code is des c ribed in detail by iSaiina et al.l 
(|2006() and iHiner et~aTl ([20091) . Here we briefly outline 
the modeling procedure. The SEDs of LoBALs are fit by 
four components: 

(1) a QSO com ponent, constructed f rom the line- free 
continuum of the iRichards et al.1 (|2006| ) SED composite 
and the emission lines from the SDSS quasar composite of 
iVanden Berk et all (I2001D (this modification was neces- 
sary to reduce the ho st galaxy contribution at lo ng wave- 
lengths present in the lVanden Berk et al.1 (|2001l ) compos- 
ite); 

(2) a hot mid-IR component, modeled as a power law 
with an exponential cutoff at short and long wavelengths, 
with a turndown at ~ 20 /im; 

(3) a warm FIR component, accounting for the small 
grain dust emission spanning a wide range of tempera- 
tures is modeled as a power law with cutoffs at high and 
low frequencies; 

(4) and a cold FIR component, modeled as a modi- 
fied black body at a fixed temperature Tdust ~ 45 K. 
Constraining the dust temperature to 45 K is a conser- 
vative assumption for a typical ULIRG-level starburst 
in nuclear regions, which are probably the closest local 
analogues to starbursts in powerful quasars. 

(5) The SEDs of the control sample require an ad- 
ditional near-infrared (NIR) component accounting for 
the emission from very hot dust, which was modeled as 
a modified black body at a temperature of T=1000 K. 
This very hot dust component, which is thought to be 
emission from the inner wall of the dust torus, is not 
necessary for fitting the SEDs of the LoBALs, with the 
possible exception of SDSS J105102+525049 and SDSS 
J170010+395554. Since we do not have NIR data for the 
majority of the LoBALs of reliable quality (see section 
§3.1.3)1 . most of their SEDs are not well constrained in 
that region and we cannot say for certain that the 3 fim 
bump is not present in LoBALs. 

The model also incl udes the Small Magellanic Cloud 
(SMC ) extinction law (jPrevot et al.lll984t iBouchet etaLl 
119851 ) and a composite PAH template or 9.7 /im silicate 
emission, if applicable. 

The FIR emission from cold and warm dust is at- 
tributed to star formation, while the mid-IR hot dust 
is assumed to be heated by the AGN. For the objects 
not detected in the FIR. 3a fluxes were used to impose 
upper limits to the SEDs. 

In Fig. [3] we show the SEDs of the 15 LoBALs for 
which both IRS spectra and MIPS photometry could be 
obtained. The SEDs of the control sample are plotted 
in Fig. 21 Overlaid are the best-fit SED model and the 
phenomenological break-up of the different components. 
Although we have observed 20 of the 22 LoBALs with 
the Spitzer IRS (see Fig. [TJ, five of those could not be 
observed with MIPS due to the early cryogen depletion. 

4.3. Measured quantities 

Our phenomenological SED modeling is used to break 
down the AGN and starburst contributions to the FIR 



emission. In order to quantify the infrared luminosities 
and the level of star formation activity in LoBALs and 
compare it to that in non-LoBALs, we estimate the fol- 
lowing quantities, which we list in Table [5] and in Ta- 
ble [6] for the samples of LoBALs and non-LoBALs, re- 
spectively. In these tables, columns (1) and (2) list the 
object number in the sample and its name. 

Column(3) gives the total infrared luminosity from 8 
to 1000 /im, i/fl^—iooo^im' which is integrated from the 
best-fit SED model and includes the contribution of the 
starburst and the AGN to the infrared flux. 

Column (4) lists the FIR luminosity contributed by 
the starburst, Lj^j. It is estimated by summing the 
warm and cold components of the model. In our SED 
model, those two components account for the emission 
from warm, small grain dust and cold dust at Tdust ~ 45 
K, respectively. We tested the robustness of the choice of 
a 45 K dust by allowing Tdust to vary as a free parameter 
in the fit, and found that it introduces a variation in the 
estimate of LfJ? of less than 5%. Attributing the warm 
and cold F IR emission to st ar formation is empirically 
justified bv lLacv et al.l (|2007| ) who find, for a sample for 
six type-2 AGN, that the sum of these components scales 
with PAH luminosity for a wide range of FIR luminosi- 
ties. Additional evidence for t he star format i on ori gin of 
the FIR emission comes from iNetzer et al.l ([20071 ) who 
find that to be the case for PG QSOs. 

Column (5) lists the contribution of the AGN (hot dust 



continuum) to the FIR luminosity, Lf^g_ 



1000pm ' 



It 



was estimated by integrating the reddened hot dust com- 
ponent and the QSO continuum from 8 to 1000 /im. 

Column (6) shows the relative percentage contribution 
of the starburst infrared luminosity to the total IR lumi- 



nosity from 8 to 1000 /mi, 



T total 
/-R,8-1000m^ 



x 100%, obtained 



as the ratio of the integrated QSO composite and hot 
dust from 8 to 1000 //m to the total integrated flux from 
8 to 1000 /im. 

In Column (7) we list the total infrared luminosity from 
3—1000 /im, Ljlj£ s . While we note that this quantity 
is calibrated for galaxies, we estimate it solely for com- 
parison with other studies. It is calculated from its em- 
pirically calibrated relationship to the bro ad-band MIPS 
fluxes at 24, 70, and 160 //m derived by iDale fc Heloul 
(l200l eq. 4): 

l ¥/r S = Ci^(24/im) + C 2 ^(70/xm)+C3^ I /(160/um) 

where [CijC2,C3] are redshift dependent coefficients (e.g., 
[3.91499, 0.48179, 1.0049] for z = 0.55). The authors 
state that this relationship reproduces the model bolo- 
metric infrared luminosities from 3 to 1000 /im to better 
than 1% for galaxies at z = 0, less than 4% for all other 
rcdshifts z < 4, and 7% for colder galaxies. The data 
used to constrain the Dale & Helou SED models ranges 
in Ltir from less than 10 s L Q to 10 12 L Q , hence, this 
method should be successful in estimating Ijtir for our 
sample of ULIRGs and LIRGs. Five of the LoBALs in 
our sample were not observed with MIPS. For three of 
those we used Scanpi 2 6 to extract IRAS photometry at 

6 Scanpi 2 is the scan processing and integration tool for ex- 
traction of IRAS photometry, developed at IPAC /Infrared Science 
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Fig. 3. — Model fits to the SEDs of the 15 LoBALs in the sample with available Spitzer IRS and MIPS observations. The solid black 
lines are the Spitzer IRS spectra an d the open black circles are SDSS ugriz and Spitzer MIPS photometry at 24, 70, and 160 fim. The 
SED model is described in detail in 94.21 The overall fit to the SED is plotted with a solid red line. The individual components to the 
fit are plotted with dot-dash lines, color-coded as follows: unreddened SDSS quasar composite in green; near-IR modified black body at 
temperature T=1,000 K in cyan; modified mid-IR power-law component in orange; warm small grain dust in magenta; and 45K modified 
black body component in blue. Down arrows indicate 3<r upper limits. Partial object names are indicated in the upper right corner of each 
plot. 
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Fig. 3. — Continued. 
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Fig. 4. — SED models of the control sample of 35 non-LoBALs. The solid black lines are the Spitzer IRS spectra and the open black 
circles are SDSS ugriz. Spitzer IRAC 3.6, 4.5, 5.8, 8.0 ^m, and MIPS photometry at 24, 70, and 160 ^m. The overall fit to the SED 
is plotted with a solid red line. The individual components to the fit are plotted with dot-dash lines, color-coded as follows: unreddcncd 
SDSS quasar composite in green; near-IR modified black body at temperature T=1,000 K in cyan; modified mid-IR power-law component 
in orange; warm small grain dust in magenta; and 45K modified black body component in blue. Down arrows indicate 3cr upper limits. 
Partial object names are indicated in the upper right corner of each plot, with the exception of the additional very hot dust component, 
plotted as dot-dashed cyan line. 
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Fig. 4. — Continued. 
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25, 60, and 100 /xm from which we estimate t heir to- 
tal IR luminosity, following iDale fc Heloul ([2002L eq. 5). 
The upper limits of Ltir are calculated using the detec- 
tions and/or the 3er flux values for the bands in which the 
source was not detected. We note that the IRAS fluxes 
may not be reliable due to the large beam size and pos- 
sible contamination by neighboring infrared sources. For 
instance, we measure the IRAS flux at 25 /jm of SDSS 
J161425+375210 with Scanpi 2 to be 130 mJy, while its 
MIPS flux at 24 (im is 20 mJy. 

Column (8) gives the SFRs calculated from the star- 
burst infrared luminosities li sted in columns (4) . The 
SFRs are estimated using the lKennicuttl (|1998l) relation- 
ship: 

SFRiMQyr- 1 ) = 4.5 x IQ- AA L IR 

where Lir is total infrared luminosity in erg/s, within 
the range 8—10 00 /urn. 

Note that the lKennicuttl (|1998l ) SFR relationship uses a 
slightly different definition of I^tir from the one derived 
with the Dale & Helou (2002) formulae, 8-1000 ^m as 
opposed to 3—1000 fj,m. The starburst components of 
the SEDs do not extend blueward of 7 /im in any object, 
so this does not affect our SFRfjf . 

For the five objects with IRS spectra that do not have 
MIPS photometry, we fit the SDSS photometry and the 
IRS data with an AGN composite and mid-IR hot dust 
component only, and, thus, estimate only the AGN con- 
tribution to the FIR. Without FIR data, it is impossible 
to estimate the starburst contribution to the FIR flux 
in those cases. However, modeling allows us to estimate 
lower limits to ^/j^g-iooo/m f rom the AGN flux from 
8—1000 /im. For all other quantities, upper limits indi- 
cate that 3cr MIPS fluxes were used to constrain the SED 
model in the FIR. 

4.4. Infrared luminosities of LoBALs 
4.4.1. Total IR luminosities from MIPS photometry 

First we address the total infrared luminosities of LoB- 
ALs to see if thei r apparent association wit h ULIRGs 
found at z < 0.4 ([Canalizo k, Stockton! 12002) is typical 
for all LoBALs. We use the definition of a ULIRG by 
iSanders et all (jl988l) . i.e., L IR > 10 12 L Q . 

For fair comparison with other studies which address 
the total infrared luminosities based on photometry mea- 
surements alone, we use the e stimates of Lj r from M IPS, 
l ¥/r S i calculated with the IDale fc Heloul (120021 ) rela- 
tion (Column (7) in Tablc[5|). Considering only the detec- 
tions, we find that one (5%) LoBAL is a LIRG (L^/£ s = 
lO n - 12 i ) and that 7/20 (35%) have L™/£ s > 1O 12 L . 
Three of those seven ULIRGs are potential HyLIRGs 
(Ljr > 10 13 L©) based on IRAS photometry, which is 
unreliable due to the large beam size and possible inclu- 
sion of neighboring sources. However, we note that more 
than half of the objects were not detected. 

Only five LoBALs in the sample have detections in the 
FIR MIPS bands at 70 and 160 /im. Of those five, only 
four are ULIRGs. Regardless of the prevalence of upper 
limits, our results unambiguously show that low-redshift 
LoBALs are not exclusively associated with ULIRGs, 
with at least 40% of them being found in LIRGs. 

Similar estimates of Lir arc derived by integrating 
the flux from 8 to 1000 /im of the best-fit SED mod- 



els (Fig. 13]). The values of i/j^g-iooo/im are listed m 
Table column (3). We find that 4/20 (20%) of the 
LoBALs are ULIRGs and 16/20 (80%) are LIRGs. For 
comparison, we estimate the total infrared luminosities, 
i/'fl's-ioooMm' °f the control sample of 35 non-LoBALs in 
the same way (see Table |6|). Of the 17 control QSOs with 
FIR detections, nine (at least 26% of the total sample) 
are ULRIGs. We find that the fraction of LoBALs re- 
siding in ULIRGs is similar to that of non-LoBALs. The 
median total infrared luminosities, i/^Liooopm' 01 the 
FIR-dctcctcd subsamples are ~ 12.12 for the LoBALs 
and 11.92 for the non-LoBALs. 

In order to compare the LoBAL and non-LoBAL sam- 
ples in the presence of so many u pper limits, we use 
the su rvival analysis statistical tests (Feigelso n fc NelsonI 
119851) . We test the hypothesis that LoBALs and non- 
LoBALs have the same distributions of infrared lumi- 
nosities. Using the Gehan and the logrank tests, we find 
that there is about 59% and 24% chance of observing 
the difference in medians if the LoBAL and the non- 
LoBAL samples were drawn from the same distribution. 
We conclude that LoBALs do not show statistically sig- 
nificant differences in infrared luminosity compared to 
non-LoBALs. 

We also find that LoBALs do not harbor intrinsically 
more infrared luminous AGN. In Fig.[5l we show the dis- 
tributions of the AGN infrared luminosities for LoBALs 
and non-LoBALs. The medians for the two types of ob- 
jects are 11.61 and 11.48, respectively. This implies that 
the infrared-to-optical ratio in LoBALs and non-LoBALs 
is comparable. This suggests that the dusty material ob- 
scuring the nuclear source has similar covering fractions. 
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Fig. 5. — LoBALs and non-LoBALs: Distribution of the AGN 
infrared luminosity integrated from the best-fit SED model between 
8— 1000 (im, plotted in bins of 0.25 dex. The medians of 11.61 for 
the LoBALs and 11.48 for the non-LoBALs show that LoBALs do 
not harbor intrinsically more infrared-luminous AGN than those 
in non-LoBALs. 



4.4.2. IR luminosity due to Star Formation 

Although in normal galaxies the FIR emission is 
usually attributed to dust emission excited by star- 
formation, it is now known that in AGN the central nu- 
clear sou rce contributes significantly to the du st heat- 
ing (e.g.. lEvans et alJl2006t1Veilleux et al.ll2009t) . Hence, 
removing the AGN contribution to the IR emission is 
important when quantifying the star formation activity 
in the galaxy and deriving star-for mation rate s from the 
total infrared luminosity using the iKennicuttl (|1998l ) re- 
lation. 

We estimate the star formation contribution to the far- 
infrared luminosity, Lj^ , by integrating the warm and 
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TABLE 5 

LoBALs: Infrared luminosities and SFRs. 



# SDSS Object ID Log(L?$_ m ) Log(Lfg) Log(L?°»_ m ) „f™ Log(L^/P) SFRfg 

[Log(L e )] [Log(L Q )\ [Log(L e )] [%] [Log(L Q )] [M yr" 1 ] 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


J023102.49-083141.2 










< 12.08 




2 


J023153.63-093333.5 










< 11.63 




3 


J025026. 66+000903. 4 


12.39 


12.25 


11.81 


73 


12.49 


310 + on 
J w —30 




T083525 984-43521 1 2 


> 11.61 








13.47 * 




5 


J085053. 12+445122.5 


< 11.87 


< 11.46 


11.65 


38 


< 12.01 


< 49 


g 


1085215 664-492040 8 


< 11.64 


< 11.52 


11.03 


75 


< 11.65 


< 56 


7 


J085357.87+463350.6 


< 11.64 


< 11.45 


11.20 


63 


< 11.69 


< 48 


8 


J101151.95+542942.7 


12.01 


11.93 


11.23 


83 


12.04 


148 +1* 


9 


J102802.32+592906.6 


> 11.04 












10 


J105102. 77+525049. 8 


< 11.72 


< 11.49 


11.32 


59 


< 11.83 


< 53 


11 


J105404.73+042939.3 


> 11.28 












12 


J112822.41+482309.9 


> 11.98 








13.39 * 




13 


J114043.62+532439.0 


> 11.77 








13.47 * 




14 


J130952.89+011950.6 


< 12.03 


< 11.56 


11.85 


33 


< 12.24 


< 62 


15 


J140025. 53-012957.0 


< 11.86 


< 11.76 


11.17 


79 


< 11.85 


< 98 


16 


J141946.36+463424.3 


11.71 


11.56 


11.18 


70 


11.84 


63 


17 


J142649.24+032517.7 


< 11.85 


< 11.49 


11.60 


13 


< 11.95 


< 52 


18 


J142927.28+523849.5 


< 12.07 


< 11.48 


11.94 


25 


< 12.27 


< 52 


19 


J161425. 17+375210. 7 


12.52 


12.28 


12.15 


57 


12.56 


326 + jg 


20 


J170010.83+395545.8 


12.12 


11.97 


11.59 


70 


12.06 


161 i?° 


21 


J170341.82+383944.7 


< 11.90 


< 11.42 


11.72 


33 


< 12.15 


< 45 


22 


J204333.20-001104.2 


< 11.84 


< 11.40 


11.64 


36 


< 11.92 


< 43 



Notes: Bold face denotes a ULIRG, i.e., L/jj > 10 Lq 

* Indicates Lj-jr estimated from IRAS rather than MIPS fluxes. 



cold components of the best-fit SED model, thus, remov- 
ing the AGN contribution. The median star formation 
luminosity of our LoBAL sample is log(Lj^ / Lq) ps 11.52 
and ss 11.35 for the control sample of non-LoBALs, when 
considering all objects, including those with MIPS up- 
per limits. LoBALs have higher median star formation 
luminosities, log(Lj^ LoBAL / Lq) w 11.97, than the non- 

LoBALs, log(Lf£ con troi/ L&) ~ 11.72, when we compare 
only the subset of objects with FIR detections for which 
the SEDs are well constrained. 

The median contribution from star formation to the 
total IR luminosity from 8—1000 /im, tota , IR (Ta- 

^IR,8 — lOOO^m 

ble[51 column (6)), in LoBALs is 41% and in non-LoBALs 
48% (70% and 62%, respectively, if only the FIR detec- 
tions arc considered) . Although we do observe significant 
variations among individual objects, from as low as 25% 
to as high as 83% in LoBALs and between 3% and 84% 
in non-LoBALs, there is no significant difference between 
the samples. In the presence of mostly upper limits in the 
FIR MIPS bands, we use the Gehan and the logrank sur- 
vival analysis statistical tests (jFeigelson fc Nelson! 1 19851) 
and find the the probability of observing the slight differ- 
ence in the median values is 79% and 96%, respectively, 
if LoBALs and the control non-LoBALs were drawn from 
the sample parent population. 

Our estimate of the star formation activity in LoBALs 
is much higher than th at found for all types of BALs by 
IGallagher et al.l (|2007l) . They report that less than 20% 
of the total FIR flux i n BAL s arises due to star forma- 
tion. Gallagh er et al.l (|2007f ) model the radio to x-ray 



SEDs of a large sample of 38 BALs at z > 1.4, consisting 
of 32 HiBALs and 6 LoBALs. Although they state that 
the quasar likely dominates the far-infrared emission in 
BALs, they note that the two most luminous starbursts 
in their sample are LoBALs, and the preponderance of 
upper limits at far-infrared wavelengths for the majority 
of the Lo- and Hi-BALs hampers their ability to draw 
definitive conclusions on the issue. 



4.5. Trends with IR luminosity 

In Fig. [6] we plot the IR luminosity from the starburst, 
Lfft, versus the AGN contribution to the IR flux from 
8-1000 /im, L^j^^_ 1 QQQ Um . At a first glance, we sec 
that the presence of silicate emission in LoBALs is cor- 
related with the AGN IR luminosity, which may suggest 
that the weak silicate emission we see in LoBALs is di- 
rectly excited by the central nuclear source. All LoBALs 
with log(Lfftg_ WQOnm / 'Lq) > 11.55 show silicate emis- 
sion. However, whether the silicate feature is detected 
is entirely dependent on the signal-to- noise ratio (S/N) 
of the IRS spectra. We tested this by introducing arti- 
ficial noise to the spectra of the objects with detected 
silicate emission using the I RAF task mknoise, creating 
spectra with S/N equal to the median of the silicate non- 
detections, S/N ~ 1.5. The silicate emission we observe 
in nine of the LoBALs would not have been detected if 
the IRS data had lower quality. We conclude that the 
non-detection of silicate in the majority of the LoBALs 
might be simply a low S/N effect. 

PAH emission is favored among the LoBALs with star- 
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TABLE 6 

Control sample of non-LoBALs: Infrared luminosities and SFRs. 



# 


SDSS Object ID 


Log(L IR u '"^ ) 




- L °S(' L /fl,8-1000/jm) 


rSB 

I R 


Log(L^/P s ) 


SFRff* 


r» — 10UU/tm 












Zi?. 

[%] 




[Mq yr" 1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


110393ft 99-1-^80033 Q 
j iuo zou . z ijouuoo . y 


< 11.47 


< 11.20 


11.14 


53 


< 11.69 


<r 97 

\ Z ( 


2 


T103333 0-94-^89818 8 


< 11 12 


< 11.00 


10.49 


76 


< 11.39 


<r 1 7 
\ 1 1 


3 


T103fi51 0-44-^7^0.^0 0. 


11.72 


11.55 


11.22 


68 


11.86 


61 


4 


T103791 1^4-^007^ 7 
j iuo ( zi.iOTjyu / oo. i 


< 11 44 


<^ 1 1 08 

<^ 11. UO 


11 19 


l 3 
io 


x 11 7^1 
<v, 11. Id 


on 

v. ZU 


5 


11 04910 9^4-^049^3 ^ 

d lUrtZlU. iJi^wf"^" 1 -' ■ O 


< 11.39 


< 11.10 


11.08 


50 


< 11.68 


- 91 

v. zi 


5 


1104^9*^ 734-^0^499 ft 


11.81 


11.50 


11.51 


19 


12.08 


54 


7 


T1 flzL^^R 844-^70747 
j lu^iouu.o'in^o iuiii . u 


< 11.32 


< 10.97 


11.06 


44 


< 11.69 


\ 1U 


Q 
O 


7104(^9^ 094-^84830 1 


\ 11. OU 


1 54 


11 44 


11 


/ 11 Q7 


U 


9 


T1fl4fi^3 704-^71 ^30 4 


< 11.67 


< 11.25 


11.47 


37 


< 11.85 


< 30 


10 


11 04840 984-563635 ft 


< 11.77 


< 11.37 


11.55 


39 


< 12.00 


< 40 


1 1 


T1 048^7 094-^01 19 3 


12.01 


1 1 03 

ii.yo 


11 20 


o 1 


1 9 07 
1Z.U I 


148 


12 


tio^ooo 914-^81004 9 


<^ iz.uy 


<^ 1 1 3ft 

<^ 11. ou 


1 9 00 
1Z. uu 


"1 8 
1 o 


<^ 1Z.OU 


v. oy 


13 


TIOMOft 194-^01(^9^ 1 
j iuoiuu. i iTJyiUiO . i 


12.19 


11.96 


11.81 


58 


12.38 


158 


14 


1105518 084-570493 5 

d IUOOIO-UOTiJ ' U4ZO.O 


< 11.80 


< 11.22 


11.66 


26 


< 12.06 


< 28 


"1 F ! 
l-> 


T105ft04 004-581 593 4 

d IUOUUI.UUTwOIOZO.'I 


1 1 Q8 

<^ ii. yo 


11. zu 


11 CO 

ii.oy 


1.7 


<^ 1 9 31 
<^ IZ.Ol 


<r^ 31 
^- Ol 


11) 


T105Q5Q 034-574848 1 

d luoyoy.yoTJ (4010.1 


11 72 


11 23 


1155 

1 l.OO 


'^9 

OZ 


1 1 Q7 
ii.y i 


90 


17 


T1 ftOOl 5 R84-559950 Q 
d iuuu i o. uoTOOzzjy . y 


11.99 


11.70 


11.68 


54 


12.16 


86 


18 


1160198 544-544521 3 


12.54 


12.33 


12.13 


64 


12.67 


367 


iy 


T1ft0341 444-541501 5 


1 1 37 


«?-* 1 1 1 8 
<v. 11. io 


in q 1 ? 
iu.yo 


uo 


<^ 1 1 

<*. 11. uo 


<^ ZO 


zu 


T1ft0593 104-545ft13 3 
d IUuOZO.IUt^OUIO.O 


1 1 ftQ 

ii.uy 


11 18 


1 1 53 

1 l.OO 


"30 
ou 


1 1 Q7 

1 1 . y / 


9 ft 
ZU 


21 


11ft0ft30 (S04-549007 5 

d IUUUOU.UUtO^^UU 1 . o 


< 12.11 


< 11.35 


12.02 


17 


< 12.41 


< 38 


22 


T1 ftOQ08 054-5331 53 9 
d luuyuo.yoTJOoiiJO . z 


< 11.92 


< 11.50 


11.70 


38 


< 12.18 


< 55 


9 3 
ZO 


T1ft3031 4ft4-410145ft 

d lUOUOl.'lUT^iUllO.U 


11 45 


11 27 


1 Q8 

iu.yo 


uu 


I 1 ftO 

I I .uu 


oz 


21 


T1 ft31 35 4ft4-40575fi 4 

d IUOIOO.'IUt^UO i OU.4: 


1 1 QQ 
J. J. .tit? 


1 1 fi5 
11. uu 


11 72 


1U 


12 22 


77 


25 


I1ft3143 7ft4-404735 ft 

d lU01'-±0. 1 UT^U4 ' oo.u 


11.65 


11.33 


11.36 


4g 


11.83 


37 


26 


J163352. 34+402115. 5 


11.96 


11.75 


11.53 


62 


12.13 


98 


27 


J163502.80+412952.9 


< 11.40 


< 10.72 


11.29 


20 


< 11.64 


< 9 


28 


J163854.62+415419.5 


< 11.70 


< 11.12 


11.56 


26 


< 11.98 


< 23 


29 


J171126.94+585544.2 


11.85 


11.68 


11.34 


68 


12.10 


82 


30 


J171334.02+595028.3 


< 11.69 


< 10.29 


11.67 


3 


< 12.11 


< 3 


31 


J171736.90+593011.4 


11.93 


11.61 


11.64 


18 


12.17 


71 


82 


J171748.43+594820.6 


11.97 


11.62 


11.71 


11 


12.12 


72 


33 


J171818.14+584905.2 


11.96 


11.80 


11.45 


69 


12.10 


109 


34 


J172104.75+592451.4 


12.06 


11.92 


11.51 


71 


12.33 


144 


35 


J172414.05+593644.0 


< 11.65 


< 11.16 


11.48 


32 


< 11.93 


< 25 



Notes: Bold face denotes a ULIRG. Objects with Lm > IO 11 ' 95 !/© were considered. 



burst IR luminosities log(Lj^ '/ 'L©) > 11.75. Higher 
total IR luminos ities for objects with P AH detections 
are also seen by ISchweitzer etHI ((2006) in PG QSOs. 
Four of the five LoBALs with MIPS detections at 70 
and 160 fj,m are ULIRGs (i.e., i^_ 1000(ttm > 1O 12 L ). 
Three of those objects show the strongest PAH fea- 
tures: SDSS J101151+542942, SDSS J161425+375210, 
and SDSS J170010+395545. There is one exception: 
SDSS J025026+000903 is a ULIRG and shows no PAHs, 
but has the strongest silicate emission feature of the en- 
tire sample. 

4.6. Star formation rates in LoBALs 

It has been shown that the total IR luminos- 
ity is a plausib le star formation r ate indicator (e.g. , 
Kennicutd[l998la f B - 8-1000 ^m): IKewlev et alj|200a 
Mann et al.l 120021 (L IR - 3-1000 fim)). However, in 
AGN, much of the accretion disk UV and optical con- 
tinuum is reprocessed by dust near the active nucleus 
and re-radiated in the infrared. By breaking down the 
different components that contribute to the overall SED 
(i.e., AGN, hot, warm, and cold dust), we partially al- 



leviate the problem of AGN contamination to the FIR 
flux. With the explicit assumption that the starburst 
component of the IR luminosity is dominated by warm 
and cold dust re-processed light from O and B stars, 
we calculate SFRs from t he total st a rburs t infrared lu- 
minosity, Lj^, using t he iKennicuttl (11998D relationship 
(Table El column (8)). iRieke et al.l (|2009D show that the 
choice of initial mass function (IMF) is crucial in cali- 
brating the SFR and, by adopting a Salpeter-like slope 
with more sha llow slope a t low masses, estimate a cor- 
rection to the IKennicutt] (|1998D SFR of SFR ffiefce09 = 
0.66 x SFRKennicutt98- Although suc h an IMF fits bet - 
ter extragalactic star forming regions (IRieke et al.ll2"00ll 
the SF Rs given here are estimated with the IKennicuttl 
(|1998|) relationship and are not to be interpreted liter- 
ally but comparatively. 

Even removing the AGN contribution to the FIR flux, 
there is possibly still contribution from older stellar 
popu lations (e.g.. lDevereux fc Ealeslll989HPopescu et al.l 
|2000() and from the AGN itself. |ffiner_elaL| (j2009D find 
evidence for non-starburst contribution to what we call 
the starburst IR flux. They model the SEDs of a sam- 
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Fig. 6. — Starburst luminosity vs. AGN contribution to the FIR 
luminosity in LoBALs and in non-LoBAL QSOs, as estimated from 
the SED models. The data are listed in Table \E\ and \E\ LoBALs 
are plotted with red circles; non-LoBALs are plotted with green 
squares. Black X denotes a LoBAL with detectable 10 fivci silicate 
emission. Filled red circles indicate the present of PAH emission, 
while the red open circles are LoBALs without PAH features. The 
presence, or absence, of silicate emission in LoBALs is correlated 
with the total AGN FIR luminosity, suggesting that the weak sil- 
icate emission we see in LoBALs is directly excited by the central 
AGN. All LoBALs with Log(Lfg^_ 1000/lrn /L e ) > 11.55, show 
silicate emission. Although the more IR-luminous AGN appear to 
preferentially have silicate emission, this is likely an artifact of the 
systematically higher S/N of their IRS spectra. PAH emission is fa- 
vored among LoBALs with starburst luminosities Log(Lj^ / Lq) > 
11.75. One exception is SDSS J025053+000903, which has the 
strongest silicate emission feature in our sample. 



pie of type-1 and type-2 QSOs in the same way we do, 
and they estimate SFRs with the IKennicuttl (|1998[) re- 
lationship using the total model IR luminosity corrected 
for the AGN contribution and the total integrated PAH 
luminosity. They find slightly lower SFRs derived from 
the PAH luminosity than from the starburst IR luminos- 
ity. This discrepancy is interpreted as the presence of an 
additional AGN contribution to the FIR flux, which has 
not been accounted for in the modeling. However, we also 
note that PAH emission can be affected by the presence 
of dust because the silicate opacity curve peaks in close 
proxi mity to some of the PAH features ([Kemper et al.l 
120041 ). And although mounting evidence supports the 
prediction that PAH carr iers are destroyed by the AGN 
radiation (e.g.. lVoitlll992D . it is still not clear whether or 
not the presence of an AGN enhances the PAH emission 
because the AGN emits UV radiat ion, which, in prin ci- 
ple, can excite PAH emission (e.g., ISmith et alJl200l . 

We find that the host galaxies of LoBALs have a range 
of star formation rates. With the caveat that most of 
our results are upper limits due to non-detections in the 
MIPS 70 and 160 /im bands, we find that the median 
SFR in LoBALs is on the order of — 52 q yr -1 and in 
non-LoBALs ~ 38 © yr" 1 . Four LoBALs have partic- 
ularly high star formation rates of ~ 150 M Q yr _1 , for 



SDSS J101151+542942 and SDSS J170010+395545, and 
-300 M yr" 1 , for SDSS J161425+375210 and SDSS 
J025026+000903. Similar fraction of the non-LoBALs 
have SFRs -100-300 M Q yr" 1 . Wc note that two of 
those four LoBALs (SDSS J161425+375210 and SDSS 
J170010+395545) and all, but one, of the highly star- 
forming non-LoBALs, do not have strong constrains on 
the FIR emission from cold dust since we only have upper 
limits at 160 /im. If we consider only the five LoBALs 
and 17 non-LoBALs with well constrained far- infrared 
luminosities, we find higher median star formation rate 
in LoBALs, -160 M Q yr" 1 , than in non-LoBALs, - 90 
M yr" 1 . 

5. DISCUSSION 
MID-IR SPECTRAL PROPERTIES 

The low-resolution mid-IR spectra of LoBALs show a 
wide range of properties. We find that about 45% of 
the LoBALs show weak 9.7 /an silicate emission (S9.7 
= 0.34—0.81), which is typical of other type-1 QSOs. 
Weak PAH emission is observed in 25% of the LoBALs, 
testifying to the presence of current star-formation in 
their hosts. Although about 40% of the LoBALs have 
featureless, low S/N spectra, we determined that the low 
quality of the mid-infrared data affects the detection of 
silicate. 

The 9.7 /im silicate feature is exclusively seen in weak 
emission in about half of the LoBALs. This supports 
the previously observed dichotomy between type-1 and 
type-2 QSOs, that is, the silicate feature appears in weak 
emission in the former and in weak absorption in the lat- 
ter. Even the ULIRGs among LoBALs exhibit silicate 
emission rather than the typical deep silicate absorption 
seen in the majority of UL IRGs (e.g.. lSpoon et al1l2007t 
sec 1 also lHao et al.l l2005a[) . Hence, in terms of their sil- 
icate feature, LoBALs resemble type-1 QSOs with their 
typical weak silicate emission. In the context of the ori- 
entation model in which BALs are present in all type-1 
QSOs, but observed only at limited viewing angles, the 
similarity of the silicate dust emission in LoBALs and in 
non-LoBALs indicates that the viewing angle for LoB- 
ALs is, not surprisingly, closer to those of type-1 QSOs 
than to those of type-2 QSOs. On the other hand, in the 
context of an evolution model in which deeply embedded 
AGN evolve to become unobscured type-1 QSOs, our 
results would imply that LoBALs mark one of the last 
stages of the transition. At this stage, the nuclear region 
has been cleared up of the thick dust envelope responsi- 
ble for the deep absorption seen in other dust-obscured 
objects such as ULIRGs. 

Yet, the optical spectra of our LoBALs suffer from high 
levels of obscuration. Wc estimate the AGN extinction 
in the optical from the SED fit, which includes SMC 
extinction law, and find a median value of Ay = 0.43 and 
a range Ay = — 1 (Table [7]), indicating significant levels 
of obscuration in some LoBALs. LoBALs are known to 
have intrin sically bluer optical-U V continua than normal 
QSOs (i.e.. lRichards et al"]l2002[ ). and so the derived A^ 
values are likely underestimated. The non-LoBALs of 
the control sample suffer significantly less extinction at 
a median level of Ay = 0.06 (Table [7]). Wc estimate 
the median color excess, E{B — V) in LoBALs, using 
the SMC extinction law and Ry = 2.72, to be E(B - 
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TABLE 7 

SMC-TYPE AGN EXTINCTION IN V BAND AND SILICATE 
STRENGTHS. 



# 


SDSS Object ID 


A(V) 


Sig.7 


LoBALs 


1 


.1023102.49-083141.2 






2 


J023153. 63-093333. 5 






3 


J025026. 66+000903. 4 


1.03 


0.81 


4 


J083525. 98+435211. 2 


0.00 




5 


J085053. 12+445122.5 


0.14 


0.34 


6 


J085215. 66+492040. 8 


0.41 




7 


J085357. 87+463350. 6 


0.53 




8 


J101151. 95+542942. 7 


0.52 




9 


,1102802.32+592906.6 


0.00 




10 


1105109 77-1-595049 8 


50 




11 


J105404. 73+042939. 3 


0.51 




12 


,1112822.41+482309.9 


0.61 




13 


Jl 14043. 62+532439.0 


0.00 
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35 
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V) ~ 0.14 a value comparable to p revi ous studies of 
LoBA Ls bv lSpravberrv fc Foltzl (|1992l ) and lGibson et al.l 
([2009h . who find SMC-type color excess of 0.12 and 0.14, 
respectively, but somewhat h igher t han the value of 0.077 
reported bv lReichard et al.1 (|2003bf ). 



INFRARED LUMINOSITIES 

LoBALs span a range of infrared luminosities. Nev- 
ertheless, they have median total and starburst infrared 
luminosities comparable to those of non-LoBALs. The 
majority of the objects were not detected in the FIR 
MIPS bands at 70 and 160 ^m. For those, roughly three- 
fourths of each sample, we have only upper limit esti- 
mates of their infrared luminosities. Using the Gehan 
and the logrank survival statistical tests, we find that 
the probabilities that the LoBAL and the non-LoBAL 
samples are drawn from the same distribution are 59% 
and 24%, respectively, considering the slight difference 
in the median total infrared luminosity, and 79% and 
96%, respectively, considering the differences in the me- 
dian starburst infrared luminosity. We conclude that the 
infrared luminosities of LoBALs are not statistically dif- 
ferent from those of non-LoBALs. 

We find a possible correlation between the extinction- 
corrected absolute Mi magnitude and the starburst lu- 
minosity of the objects with FIR detections (see Fig.[7K). 
This implies that the star formation rates are higher in 
the more optically luminous sources. However, since the 
majority of the measurements are upper limits (22 de- 
tections vs. 28 upper limits), we are not able to properly 
quantify the slope or the scatter in the correlation. Simi- 
larly, evidence for a correlation between the optical 5100 
A luminosity and the 60 IRAS flux wa s found in 
low-redshift PG QSOs bv lNetzer et al.l (|2007t ). who pro- 
posed the correlation was mostly due to star formation. 
We note that the control sample of non-LoBALs is well 
matched to the absolute Mi magnitudes of the LoBALs 
(Figure [7]) . The range of Mi, corrected for Galactic ex- 
tinction, for the LoBALs is -22.41 > Mi > -25.55, with 
a median of Mi = —24.10 (median of the FIR-detections 
is -23.68). For the control non-LoBALs, the range is 
-22.17 > Mi > -25.53, with a median of M l = -23.96 
(median of the FIR-detections is -23.85). 

We also find a correlation between the AGN infrared 
luminosity and the optical i-band luminosity (Fig. [7p), 
LfgN _ L a with a slope a ^ o.29±0.05. This relation- 
ship holds for both LoBALs and non-LoBALs and cov- 
ers over three orders of optical magnitude. If the AGN 
infrared luminosity arises mostly from dusty obscuring 
material close to the central source, then this correlation 
implies that LoBALs and non-LoBALs have comparable 
covering fractions. 

STAR FORMATION 

Star formation contributes a median of 41% and 48% 
of the FIR power in LoBALs and non-LoBALs, respec- 
tively, with large variations among in dividual objects, 
in agr eement with the value found by iSchweitzer et all 
(|2006h for PG QSOs at z < 0.3 who claim that star for- 
mation contributes at least 30% of their FIR luminosity. 

We estimate SFRs from the FIR luminosity solely con- 
tributed by the starburst. With the exception of four 
LoBALs with SFRs - 150 - 300 M Q yr" 1 , and with 
the caveat that we only have upper limits for most of 
the other objects, we find that LoBALs have SFRs ~ 52 
Mq yr _1 , slightly higher but comparable to the value we 
find for non-LoBALs, 38 M Q yr _1 . However, we note 
that for the most star-forming galaxies, those with FIR 
MIPS detections, the median SFR in LoBALs (~ 161 
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Fig. 7. — Infrared luminosity vs. absolute M; magnitude cor- 
rected for AGN reddening with the extinction estimated from the 
SED fitting. The LoBALs are plotted with open red circles; the 
non-LoBALs are shown with open green squares. Down arrows in- 
dicate upper limits. The FIR limits for the non-LoBAL QSOs are 
better constrained due to the lower uncertainties in their MIPS im- 
ages (i.e., non-LoBALs: cr7o M m ~ 1.6 mjy and eri60/jm ~ 3.7 mjy; 
LoBALs: o7o Mm ~ 6.2 mjy and crieoum ~ 8.1 mjy). (a) Starburst 
infrared luminosity vs. Mi. (b) AGN infrared luminosity from 
8—1000 (im vs. Mi. The apparent correlation implies comparable 
IR-to-optical ratios, hence, similar cover fractions in LoBALs and 
in non- LoBALs. (c) Total infrared luminosity from 8—1000 fira vs. 



Mq yr 1 ) is higher than the median SFR in non-LoB AL 
type-1 hosts (~ 90 M Q yr^ 1 ). 

LOBALS HITHERTO AND FUTURE WORK 

Overall, our results suggest that LoBALs are very sim- 
ilar to non-LoBALs in terms of their mid-infrared spec- 
tral properties and far-infrared luminosities. This fits 
the orientation scenario, which would predict compara- 
ble levels of star formation for LoBALs and non-LoBALs 
if both are drawn from the same parent population. 
Our results, however, cannot rule out an evolutionary 
paradigm, where LoBALs arc rapidly transitioning from 
a dusty phase marked by high SFRs to a more quiescent 
phase with SFRs typical of non-BALs. 

The majority of the LoBALs in our sample have SFRs 
comparable to non-LoBALs, which in the framework of 



the evolutionary model implies that most of the LoBALs 
have already passed through the event that quenched 
the star formation in the galaxy to the levels seen in 
normal QSO hosts. Our results show that there are large 
variations among individual LoBALs. Finding that only 
20% of the LoBALs have SFRs 80% higher than those 
found in the IR-luminous non-LoBAL implies that the 
period during which the star formation was quenched was 
very brief during the short LoBAL transition phase. If 
the SFRs in both samples were found to be much higher, 
but still comparable, a stronger claim could be made that 
BAL outflows alone are not responsible for quenching 
the star formation on galactic scale, but that was not 
observed. 

While the Spitzer observations explore the relation be- 
tween LoBALs and star formation, they are not sufficient 
to test whether the LoBAL phenomenon is indeed related 
to the early stages of QSO activity. Moreover, due to the 
preponderance of upper limits, our Spitzer results do not 
allow us to draw strong conclusions about the nature 
of LoBALs. This study will be complemented by HST 
imaging and Keck spectroscopy programs of this sam- 
ple, which will help us study the host galaxy morpholo- 
gies and the ages of the dominant stellar populations. 
On one hand, correlating the merger stage with spectral 
characteristics will give us a deeper insight into the dy- 
namics involved and help us constrain an evolutionary 
connection between unobscured QSOs and LoBALs. On 
the other hand, the stellar ages will help us constrain 
the time scales involved in fueling of the nucleus and the 
onset (and perhaps quenching) of star formation. 

6. CONCLUSION 

We investigate the mid- and far-IR properties of a 
volume-limited sample of 22 low-ionization broad absorp- 
tion line QSOs within the redshift range 0.5 < z < 0.6. 
We model their SEDs from the optical to the far-infrared 
in an effort to estimate total infrared luminosities, the 
relative contributions from the starburst and the AGN, 
starburst luminosities, and star formation rates corrected 
for the AGN contamination of the FIR emission. We 
compare this LoBAL sample to a control sample of non- 
LoBALs, matched by Mj within the redshift 0.45 < z < 
0.83, to examine the possible connection between these 
two classes of QSOs. We find that LoBALs are indistin- 
guishable from non-LoBAL type-1 QSO in terms of their 
MIR spectral properties and FIR luminosities. 
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APPENDIX 

APPENDIX: SDSS SPECTRA OF THE 22 LOBALS 

Figure [8] shows the SDSS spectra of the 22 LoBALs in the sample and details on the Mg II line to assure the 
reader that the sample comprises of bona fide LoBALs. We adopt the more inclusive definition bv lTrump et al.l ([2006D 
requiring the troughs to span a velocity width of at least 1000 km s _1 blue- ward of Mg II A2800. Eight of the 22 
LoBALs here have blue-shifted Mg II absorption line wi dths between 1000 and 2000 km s , which would classify them 
as mini-BALs according to the traditional definition bv lWevmann et al.l (119911 ) requiring widths greater than 2000 km 
d-1 
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Fig. 8.— Left: SDSS spectra of the 22 LoBALs. Right: Velocity profile of the Mg II A2800 line. 
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